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Abstract 


Determination of pier-scour potential is an important 
consideration in the hydraulic analysis and design of highway 
bridges that cross streams, rivers, and other waterways in the 
United States. A primary goal of ongoing research in the field 
of bridge scour is to improve scour-prediction equations so 
that pier-scour depth is neither underpredicted nor excessively 
overpredicted. Scour depth for piers in noncohesive, non- 
uniform streambeds with a mixture of sand, gravel, cobbles, 
and boulders (coarse-bed streams) generally is less than the 
scour depth in finer-grained (mostly sand) streambeds under 
similar hydraulic conditions. The difference in scour depth 
is attributed to formation of an armor layer. Pier-scour data 
collected by the U.S. Geological Survey were used to develop 
a bed-material correction factor called K,. The equation 
recommended by the Federal Highway Administration for 
computing pier scour is a version of the HEC-18 pier-scour 
equation that includes K,, which is referred to in this report 
as the HEC-18-K,Mu equation. The Montana Department 
of Transportation was interested in pier-scour prediction in 
coarse-bed streams because coarse-bed streams are common 
in Montana. Consequently, the U.S. Geological Survey and the 
Montana Department of Transportation began a cooperative 
study in 2001 to investigate pier scour in coarse-bed streams 
in Montana. 

This report describes results of a study of pier scour in 
coarse-bed streams at 59 bridge sites during 2001-07 in the 
mountain and foothill regions of Montana. Drainage areas for 
the streams at bridge sites where measurements were collected 
ranged from about 3 square miles (mi*) to almost 20,000 mi’. 
Data collected and analyzed for this study include 103 coarse- 
bed pier-scour measurements. The report also describes how 
coarse-bed pier-scour measurements were collected, shows 
the extent that the coarse portion of the national pier-scour 
database was expanded, discusses how these new data were 
used to evaluate the relative accuracy of various equations for 
predicting scour in coarse-bed streams, and demonstrates how 
differences in size and gradation between surface bed mate- 
rial and bed material underlying the surface layer (shallow- 
subsurface bed material) might relate to pier scour. 


Ninety-six out of 103 pier-scour measurements were 
made under clear-water scour conditions, when the streambed 
upstream from the bridge is stable and there is no substantial 
incoming sediment supply to the bridge opening. Of the mea- 
surements made, 50 percent had an approach velocity (V) that 
equaled or exceeded 70 percent of the critical velocity (V,.,) 
for incipient motion of bed material, which might indicate that 
scour was measured very near the threshold between clear- 
water and live-bed scour (V,/V,, equal to 1.00) where maxi- 
mum scour was shown in laboratory studies. 

Pier-scour data collected for this study were compared to 
selected pier-scour data from the Bridge Scour Data Manage- 
ment System (BSDMS), a database of scour measurements 
made nationwide, to show the effect of bed-material size 
and gradation on scour depth. The relation between relative 
pier scour (y’./b), or the ratio of measured pier scour (y’ ) to 
effective pier width (4), and relative velocity (V\/V_.,) defined 
by an envelope curve for data collected for this study and the 
BSDMS data were compared to relations observed in labora- 
tory data. The envelope curve for data used for this study 
displays an earlier peak, with a decrease in relative pier scour 
followed by an asymptotic rise. Interestingly, the first peak 
in the curve occurs at a lower relative velocity (V/V... about 
0.75) than the clear-water and live-bed threshold (V//V., 
equal to 1.00). Also, the magnitude of this peak (maximum of 
1.25) is much less than the published maximum value of 2.4 
obtained from laboratory data. Furthermore, equilibrium live- 
bed scour is about 18 percent shallower than the maximum 
pier scour under clear-water conditions, a larger reduction 
than the 10 percent indicated in published research. Armoring 
associated with nonuniform coarse-bed material and unsteady- 
flow conditions in the field probably accounts for differences 
between envelope curves developed for this study and those 
from previous studies conducted in the laboratory using fine- 
grained material and steady-flow conditions. 

Scour depth was computed for the measurements col- 
lected for this study using the HEC-18 equation without the K, 
correction factor and five pier-scour equations that use K, for 
the armoring effect of coarse bed material. The HEC-18-K,Mu 
equation was the best equation for predicting pier-scour 
depth in coarse-bed streams because the equation generally 
predicted pier scour in closer agreement to measured scour 
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than the other equations used for computing pier-scour depth. 
Furthermore, the magnitude of the residuals of underpre- 
dicted and overpredicted scour depth was the lowest using the 
HEC-18-K,Mu equation. 

Gradation coefficients (o) for the surface bed-material 
data collected for this study ranged from 1.48 to 4.14, with a 
median of 2.01. For data used in this study, maximum relative 
plier scour was greater when o, was less than or equal to about 
2.5. Pier-scour measurements collected for this study and 
BSDMS confirm a general lack of data for coarse-bed sites 
with o, greater than about 2.5. 

Paired samples of surface and shallow-subsurface bed 
material collected for this study were analyzed and the D,, and 
D,, particle sizes were compared. Median D,, particle sizes 
for 103 surface and shallow-subsurface bed-material samples 
were about 49 millimeters (mm) and 32 mm, respectively, and 
median D,, particle sizes for 103 surface and shallow-subsur- 
face bed material samples were 122 mm and 86 mm, respec- 
tively. The median gradation coefficients for surface and shal- 
low-subsurface bed material associated with measurements 
made for this study were 2.01 and 4.14, respectively. The 
frequency and magnitude of the residual pier-scour depth asso- 
ciated with underpredicted and overpredicted scour was about 
the same regardless of whether surface or shallow-subsurface 
particle-size data were used with the HEC-18-K,Mu equation. 

The combined effect that surface and shallow-subsurface 
bed-material characteristics might have on measured pier 
scour also was examined. The combined influence of relative 
differences in surface and shallow-subsurface bed-material 
size and gradation on pier scour was investigated by defining 
a bed-material variability index (/,) based on ratios of particle 
sizes and gradation coefficients between layers for each of 
the 103 pier-scour measurements. Selected percentiles were 
determined and the relative pier scour and relative velocity for 
each measurement in the percentile were plotted. Envelope 
curves bounding the measurements in each percentile were 
drawn. Maximum relative pier scour associated with each 
envelope curve was concentrated at about the same relative 
velocity. Based on these curves, as differences in particle size 
and gradation between the surface and shallow-subsurface 
bed-material layers increase (higher values of J), pier scour 
decreases. 

Runoff conditions during the study mostly were lim- 
ited to bankfull discharge approximating the 1.5- to 2-year 
recurrence-interval flood, which have a 67- and 50-percent 
chance, respectively, of happening in any year. Lack of higher 
flows having greater velocity intensities may indicate that 
scour depths were limited by hydrologic conditions. Pier-scour 
depths measured in this study also may have been limited by 
streambed armoring and unsteady flow conditions. Conclu- 
sions presented for this study are generally limited to the range 
of hydraulic conditions and bed-material characteristics dem- 
onstrated in the data collected for this study and the BSDMS 
data. 


Introduction 


Determination of pier-scour potential is an important 
consideration in the hydraulic analysis and design of highway 
bridges that cross streams, rivers, and other waterways in the 
United States. Existing bridges continue to undergo scour 
assessments as part of a national program, and pier scour is 
among the key scour components investigated. Bridge reha- 
bilitation and the design of new bridges require that pier foot- 
ings be located below the scour zone of the streambed during 
large floods. Pier-scour depth is commonly computed using 
empirical equations, with many different equations having 
been developed over the years. Such equations were tradition- 
ally based on studies conducted in hydraulic laboratories using 
flumes and scaled pier models. Scale effects associated with 
model studies restricted the maximum bed-material particle 
size to mostly noncohesive sand. A primary goal of ongoing 
bridge scour research is to improve laboratory-based scour- 
prediction equations using field-scale data so that pier-scour 
depth is neither underpredicted nor excessively overpredicted. 
Underpredicted pier scour can result in pier footings not being 
located deep enough in the streambed to resist scour. Scour 
during large floods could, therefore, undermine pier footings 
and result in damage or instability to the bridge and possibly 
catastrophic failure. Overpredicted pier scour might indicate 
that more bridges are scour critical, or prone to failure due to 
scour, than is actually the case and thus may lead to expensive 
over-design or unnecessary retrofitting of pier foundations 
(Chase and Holnbeck, 2004). Research efforts, therefore, 
continue to evaluate and improve the accuracy of pier-scour 
equations. 

Since the early 1990s, research has expanded to include 
field measurements of pier scour during high-flow conditions. 
These scour measurements enable pier-scour equations to be 
evaluated beyond the theoretical or limited laboratory condi- 
tions frequently used in equation development. Mueller (1996) 
and Mueller and Jones (1999) studied bed-material size and 
gradation effects on pier-scour depth in noncohesive stream- 
beds. Their research indicated that scour depth for piers in 
noncohesive, nonuniform streambeds with a mixture of sand, 
gravel, cobbles, and boulders (coarse-bed streams) generally is 
less than the scour depth in finer-grained (mostly sand) stream- 
beds under similar hydraulic conditions. The difference in 
scour depth is attributed to formation of an armor layer (Shen 
and others, 1966). This finding, supported by field measure- 
ments collected by the U.S. Geological Survey (USGS), has 
led to modifications in the scour equations recommended by 
the Federal Highway Administration (FHWA). Pier-scour data 
collected by the USGS were used to develop a bed-material 
correction factor called K,. Collection of additional coarse- 
bed pier-scour data was recommended to further validate the 
method for computing K, and to provide additional data for 
continued evaluation of pier-scour equations (Mueller and 
Jones, 1999). As recently as 2005, Mueller and Wagner (2005) 
noted that relations and equations developed from labora- 
tory research had not been adequately verified with field data. 


Much advancement has been made since the mid—1990s in 
the research and development of equations that account for 
coarse-bed pier scour, but additional data were needed to 
evaluate whether these equations accurately predicted scour in 
coarse-bed streams. 

The Montana Department of Transportation (MDT) 
was interested in pier-scour prediction in coarse-bed streams 
because coarse-bed streams are common in Montana. Further- 
more, the USGS has collected bridge-scour data and con- 
ducted bridge-scour investigations in Montana since the early 
1990s in cooperation with MDT. Consequently, the USGS and 
MDT began a cooperative study in 2001 to investigate pier 
scour in coarse-bed streams in Montana. 


Purpose and Scope 


This report presents the results of a study of pier scour in 
coarse-bed streams at 59 bridge sites during 2001-07 in the 
mountain and foothill regions of three major river basins in 
Montana: the Missouri, Yellowstone, and Clark Fork of the 
Columbia (fig. 1). Additionally, the report: 


1. Describes how coarse-bed pier-scour measurements 
were collected. 


2. Shows the extent to which the coarse portion of the 
national pier-scour database used for conducting 
scour research was expanded. 


3. Discusses how these new data were used to evalu- 
ate the relative accuracy of various equations for 
predicting scour in coarse-bed streams. 


4. Demonstrates how differences in size and grada- 
tion between surface bed material and bed material 
underlying the surface layer, referred to in this report 
as shallow-subsurface bed material, might relate to 
pier scour. 


Drainage areas for the streams at bridge sites where 
measurements were made ranged from about 3 mi? to almost 
20,000 mi?. Data collected and analyzed for this study include 
103 coarse-bed pier-scour measurements and represent a sub- 
stantial sample of coarse-bed streams in Montana. When these 
data are added to the USGS Bridge Scour Data Management 
System (BSDMS) database at the conclusion of the study, the 
number of coarse-bed pier-scour measurements in the current 
(2011) database will increase substantially. 

Field methods (Jarrett and Boyle, 1986; Landers and 
Mueller, 1996; and Mueller and Landers, 2000), hereinafter 
referred to as limited-detail field procedures, and equipment 
ranging from simple scour rods to complex hydroacoustic 
instrumentation were used to obtain pier-scour measure- 
ments. Pier-scour measurements were made during spring 
snowmelt runoff when streamflow and hydraulic conditions 
generally approximated bankfull discharge with recurrence 
intervals of 1.5 to 2 years, which have a 67- and 50-percent 
chance, respectively, of happening in any year. Pier-scour 
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measurements also were made during low-flow conditions. 
Because the low-flow measurements are considered ancillary, 
only the high-flow measurements are presented and discussed 
in the report. Selected Montana pier-scour measurements for 
peak streamflows (5— to about 100—year recurrence interval) 
obtained prior to this study also were used as an independent 
check on envelope curves developed similar to those from 
previous studies (Ettema, 1980; Melville, 1984; Chiew, 1984; 
and Mueller and Wagner, 2005). The envelope curves have 
either linear or curvilinear shape and define an upper boundary 
to the data. 

In this report, gradation and gradation coefficient (, ) are 
terms frequently used to describe and measure the variation in 
particle-size distribution of the bed material. Bed material is 
considered to be poorly graded or has a smaller gradation (also 
referred to as well sorted) when the material is more uniform, 
or has a narrow range of grain sizes. Bed material is consid- 
ered to be well graded or has a larger gradation (also referred 
to as poorly sorted) when the material is nonuniform, or has 
a wider range of grain sizes. The magnitude of the gradation 
coefficient increases as gradation increases or as the range in 
particle sizes increases. 

Because bed armoring relates to surface and subsurface 
particle-size differences, surface and shallow-subsurface bed- 
material data were collected at each site. Particle-size grada- 
tion was determined from pebble counts and sieve analyses. 
Surface and shallow-subsurface bed material sampled during 
this study were used to evaluate the sensitivity of computed 
scour to size and gradation differences between the two layers. 
Data for each layer also were used to assess the effects of bed- 
material characteristics on measured pier scour. 

Mathematical symbols (=, ~, <, >, <, and >) used in this 
report are defined near the front of the report under the section 
“Symbols and variables used in this report,” and are either 
spelled out or appear as the symbol thereafter in the body of 
the report. Variable symbols used in this report (for example, 
y,, D,, and V,) are defined near the front of the report under 
the section “Symbols and variables used in this report,” and at 
first occurrence and occasionally thereafter in the body of the 
report. 


Background 


Hydraulic Conditions and Real-Time Scour 
Measurements 


The relative stability of a streambed is related to hydrau- 
lic conditions in the stream and the particle sizes of the 
streambed material. Scour at a bridge pier depends on the 
stability of the streambed near the bridge, the extent of bed- 
material transport from upstream, and forces acting on the 
streambed near the pier that result from the interaction of the 
stream and the pier. Pier scour can occur when the streambed 
upstream from the bridge is stable and there is no substantial 
incoming sediment supply to the bridge opening (clear-water 
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Figure 1. Location of bridge sites with pier-scour data used for this study. 


scour) and when the streambed is unstable and bed material is 
actively being transported by the stream into the bridge open- 
ing from upstream (live-bed scour). Documenting the stream 
stability and sediment-transport conditions associated with 
each pier-scour measurement is necessary to improve under- 
standing of the processes contributing to pier scour. In order 
to develop equations relating scour and hydraulic conditions, 
pier scour measured in the field must be the result of hydraulic 
conditions at the time of measurement (real-time scour). 

The critical velocity for incipient motion of bed material 
(V.,) is used to determine if clear-water or live-bed conditions 
predominate upstream from a bridge site (Richardson and 
Davis, 1995, p. 12-13), with observed velocity less than V,, 
indicating clear-water scour and observed velocity greater than 
V_ indicating live-bed scour: 


V, =31.080"K,y,D/, (1) 


where 
0 is the Shields parameter, dimensionless, based 
on particle size and hydraulic conditions; 
K is a constant equal to 1.81 for inch-pound 
units, dimensionless; 
y, is the approach depth of flow just upstream 
from the pier, in feet for this equation; and 
D. is the bed-material particle size D, in feet for 
this equation, where i percent of particles 
are finer by weight. 
Equation | is commonly applied to 
determine the critical velocity for any 
particle size in the streambed. For 
incipient motion, the median or D,, 
particle size and corresponding 6 typi- 
cally are used to determine stream sta- 
bility. Clear-water scour often occurs 
in streams with coarse bed material. 
Pier scour is a local scour 
process resulting from vortex action 
and flow acceleration when flow is 
obstructed by the pier. The vortex 
systems produced are complex three- 
dimensional flow fields that can act 
as continuous or intermittent energy 
systems that vary with hydraulic 
conditions. The two primary vortices 
involved are the horseshoe vortex and 
the wake vortex (fig. 2). The horseshoe 
vortex forms near the upstream base of 
the pier as flow impinges and moves 
down the pier face and along the sides 
of the pier. The wake vortex forms as 
shearing layers of flow detach from 
the sides of the pier and move down- 
stream away from the pier, which 
results in bed material being carried or 
scoured away (Melville, 1995). 


Flow-acceleration 
zone 


Detailed 
perspective 


Figure 2. 
Melville, 1995). 
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The equation by Gao and others (1993) is used to com- 
pute the critical velocity to initiate scour of bed material at the 
pier (V”,, ) based on the critical velocity for incipient motion of 
bed material upstream from the bridge (V_, , equation 1): 


0.053 
y_,=0.645|/— | V,, (2) 
b, 
where 
b, is the nominal pier width measured 
perpendicular to the pier length, which 
takes on the units of D, in this equation 
(feet); 
and where all other terms are as previously defined. By this 
equation, the velocity to initiate scour of bed material at the 
pier is about one-half the velocity for incipient motion of bed 
material in the channel upstream from the bridge. 

Equation 2 is used to determine if a valid real-time mea- 
surement of pier-scour was made in the field. Measurements 
(scour depth, flow velocity, and flow depth) made at remnant 
scour holes under base or moderate flows do not reflect the 
true hydraulic conditions accounting for the geometry of the 
scour hole. The measured approach velocity (V,) is compared 
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to the critical velocity (V’,.,) to initiate scour of the D,, particle 
size computed by equation 2, to determine if the hydraulic 
conditions at the time of measurement are sufficient to cause 
pier scour, according to the following criteria: 


1. IfV>V'..,, a valid real-time scour measurement was 
made because bed material can scour at the pier under the 
observed conditions, or 


2. IfV<V"..,, a valid real-time scour measurement was not 
made because bed material cannot scour at the pier under 
the observed conditions. 


Under clear-water conditions, scour continues to develop 
until the scour hole becomes armored, and the resulting scour 
is referred to as maximum-local-pier scour. Under live-bed 
conditions, scour continues to develop until the rate of sedi- 
ment transport out of the scour hole equals the rate of sedi- 
ment transport into the scour hole, or the scour hole becomes 
armored and the stream velocity is insufficient to transport 
the size of bed material in contact with flow (Richardson and 
Davis, 2001). Scour depth under live-bed conditions tends to 
fluctuate around a mean-depth value over time, partly due to 
the transient nature of bedload transport in coarse-bed streams 
and passage of dunes in sand-bed streams. This mean-depth 
value is referred to as the equilibrium local live-bed pier scour. 
According to Richardson and Davis (1995, p. 16), the equilib- 
rium local live-bed pier scour is about 10 percent shallower 
than the maximum-local-pier scour under clear-water condi- 
tions. Maximum-local-pier scour may be attained under clear- 
water conditions before the threshold for incipient motion is 
reached in the upstream channel and before live-bed scour 
conditions result in transport of bed material into the scour 
hole. Exceptionally high flood discharges that induce high 
bed-material transport upstream from the bridge, therefore, 
may not produce the deepest scour at a pier. Local pier scour 
is frequently used to describe scour caused by local hydraulic 
conditions at the pier, where the effects of any contraction 
scour, abutment scour, general scour, or long-term scour at the 
bridge are excluded. Hereinafter, local pier scour is referred to 
as pier scour. 

Clear-water and live-bed distinctions have particular 
importance in laboratory research where either condition can 
be maintained throughout the hydraulic experiment. However, 
these distinctions become less clear under field conditions. 
During floods, piers in coarse-bed streams can be subjected 
to clear-water scour at lower flows, live-bed scour at higher 
flows, and clear-water scour again during the receding flood 
stages (Mueller, 1996). Furthermore, Mueller (1996) indicated 
that only a relatively short interval of live-bed scour is needed 
for pier scour to attain an equilibrium scour depth typical of 
live-bed conditions. Finally, Mueller and Wagner (2005, p. 19) 
noted that differences in armoring rates between the upstream 
channel and the pier-scour hole can complicate scour-hole 
development, overall scour depth, and the actual condition 
(clear-water or live-bed). They concluded that the traditional 
clear-water and live-bed designations are insufficient to fully 
describe the actual conditions for streams with nonuniform 


bed material. However, clear-water and live-bed scour des- 
ignations continue to be important concepts describing the 
predominant scour condition at a site. 


Development of Pier-Scour Equations for 
Coarse-Bed Streams 


Pier-scour research in the mid—1970s resulted in a pier- 
scour equation, which was recommended by FHWA (Richard- 
son and others, 1975) and widely adopted for conducting pier- 
scour analysis and design in the United States. This equation 
was largely based on laboratory flume studies with no verifica- 
tion using actual field data. Additional research based on field 
data in later years indicated that pier scour in nonuniform, 
coarse, bed material was not well predicted by the recom- 
mended equation. The equation has evolved to include correc- 
tion factors for pier shape, flow alignment, bed-form condi- 
tions, and bed-material size and gradation. The latest (2011) 
version of the equation incorporates a bed-material correction 
factor (K, ) to compute pier scour in coarse-bed streams. 

The pier-scour equation recommended by the FHWA for 
conducting detailed bridge-scour investigations in the mid— 
1970s was the Colorado State University (CSU) equation 
(Richardson and others, 1975). A modified version of the CSU 
equation was later published by FHWA in Hydraulic Engineer- 
ing Circular No. 18 (HEC-18) (Richardson and others, 1991) 
and is referred to in this report as the basic HEC-18 equation. 
The basic HEC-18 equation is: 


y, = 2.0y,K,K,|— | FO? (3) 
where 
y, is the computed pier-scour depth that assumes 
the units of y, (in feet for this study); 
K, is the correction factor for pier-nose shape, 
dimensionless; 
K, is the correction factor for flow angle of attack 


on the pier (a), dimensionless; and 
Fis the Froude number directly upstream from 
the pier, dimensionless; 
where the Froude number is computed according to: 


oe oe (3a) 


where 
g is the acceleration due to gravity, in feet per 
second squared for this study; 


and where all other terms are as previously defined. Further 
research produced a second version of the basic HEC-18 equa- 
tion recommended by FHWA (Richardson and others, 1993, 

p. 39) that included a correction factor (K,) to account for the 
bed-form condition (dunes, plane bed, and other variations): 


b, a 0.43 
y, =2.0y,K,K,K,) — i ‘ (4) 
where 
K, is the correction factor for bed-form 


condition, dimensionless; 


and where all other terms are as previously defined. 

The value of K, is equal to 1.1 for clear-water scour and plane- 
bed conditions—typical of coarse-bed streams analyzed in this 
study. Other values of K, range from 1.1 to 1.2 for beds with 
medium dunes to 1.3 for beds with large dunes (Richardson 
and others, 1993, p. 39). These first two versions of the basic 
HEC-18 equation did not account for bed-material size and 
gradation effects related to armoring of the streambed. 

Most of the pier-scour equations that have been devel- 
oped were based on laboratory research that used bed material 
of uniform sand or fine gravel. Consequently, these equations 
do not adequately describe the effect of bed-material size and 
gradation on scour depth under actual field conditions. McIn- 
tosh (1989) reported that more than 35 equations have been 
proposed to compute pier scour. Mueller (1996) found that of 
these 35 equations, only 8 included an explanatory variable 
(median particle size or D,,) to account for bed-material size. 

The effect of nonuniform bed material on scour depth 
was observed by Shen and others (1966), where the presence 
of larger particle sizes appeared to reduce scour depth due to 
armoring. Research by Nicollet and Ramette (1971) com- 
pared scour for three separate sand sizes and a sand mixture 
of equal parts of each size. They concluded that the maximum 
scour occurred at the critical velocity of each sand size used 
and that scour depth was less for experiments with the sand 
mixture. The bed-material gradation coefficient (9,) was used 
by Ettema (1976, 1980) to show that nonuniform bed material 
affected the rate of development and final depth of scour under 
clear-water conditions. These findings were extended to live- 
bed scour conditions by research conducted by Chiew (1984) 
and Baker (1986). Furthermore, research by Abdou (1993) 
substantiated that an increase in oO, resulted in a decrease in 
clear-water scour. In all of these studies, findings were based 
only on laboratory investigations with no validation with field 
data. 

The need to account for bed-material size and gradation 
effects on pier scour led to research sponsored by FHWA, 
where laboratory data were used by J. Sterling Jones to 
develop a correction factor (K,,) that was first presented by 
Richardson and Davis (1995) and described later by Mueller 
and Jones (1999). K, accounts for armoring of the scour hole 
and generally reduces computed scour in coarse-bed materials 
compared to what would be computed if the factor were not 
applied. 

In response to the work of Jones, the FHWA recom- 
mended a new version of the HEC-18 equation for use for 
both clear-water and live-bed scour conditions that incorpo- 
rated the K, correction factor (Richardson and Davis, 1995): 
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0.65 


0.43 
Fe, (5) 


where all terms are as previously defined. Richardson and 
Davis (1995) defined K, as: 


is 


k,=| 1-089(1-7,) "| (6) 


where 
Vi is the velocity-intensity ratio, dimensionless, 
computed as: 


V, — V ‘i 
F, 90 “a a 


c 


i= ; (7) 


where all other terms are as previously described. 


The laboratory research by Molinas and Abdou (1999) 
and Molinas and Noshi (1999) proposed the following equa- 
tions for computing K, : 


0.75 
D, V-V, Vi -V, 
K,=2.0+1.25, /—=| —-—* | in] —*>—*- |, (8) 
Ds, Vice a a Veg = Vn 


D.,, 18 the particle size that represents the coarse fraction of the 
bed material, in feet, computed as: 


Dy + 2Dy. + 2Dy5 + Dayo 
(2) cg 6 


, (9) 


V. is the critical velocity that initiates scour 
of bed material at the pier, and equal to 
0.4V_,. , in feet per second; 
Vo, is the critical velocity for incipient motion of 
the D_,, particle size, in feet per second; 
and where all other terms are as previously defined. 

Mueller (1996) used field data collected by the USGS 
throughout the United States to develop a K, factor. Partial 
residual analysis and multiple linear regression were used to 
identify the important variables affecting pier-scour depth. 

The investigation showed that the D,. particle size was the 
dominant bed-material variable affecting pier scour under both 
clear-water and live-bed scour conditions. The D,, particle 
size also was considered important based on Gao and others 
(1993), where the velocity-intensity ratio based on only the 
D,, particle size was shown to have an influence on pier-scour 
depth. Critical velocities associated with these two important 
bed-material variables resulted in the Mueller velocity-inten- 
sity ratio (V,,, dimensionless): 


Vix ~ V, = Vso , 
V 50 =i 95 


and where all other terms are as previously defined. 


M ? 


(10) 
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The following equation incorporating the variables that define 
V_.,, Was then proposed by Mueller (1996) for computing K,: 


RM 
V y' 0.15 
K=04|\— | 
r 50 V 95 


(11) 
Equations 10 and 11 apply for the condition when: 
us a Vso 
0<}| —*— |< 450. (12) 
Vso — Vios 


In a review of equations describing pier scour, Mueller 
and Jones (1999) concluded that equations 10-12 proposed by 
Mueller (1996) for quantifying K, provided the best expla- 
nation of the differences between observed and computed 
pier-scour depth. However, the authors recognized the need 
that additional field data be obtained and further research be 
conducted to improve the accuracy of pier-scour equations. 

Chase and Holnbeck (2004) compared 42 pier-scour 
measurements made by the USGS in 5 States at 11 bridge sites 
with scour depth computed by 5 of the more commonly cited 
equations. The version of the HEC-18 equation that includes 
K, as proposed by Mueller (equations 5 and 11) was statisti- 
cally determined to be the best equation for predicting pier 
scour in coarse-bed materials for that data set. 

Mueller (1996, p. 115) found that the HEC-18 equation 
(equation 5) underestimated very few instances of observed 
scour but substantially overestimated observed scour depths 
because of the restrictions placed on K, . When first published 
as part of the HEC-18 equation (Richardson and Davis, 1995, 
p. 36), K, could only be applied when the D,, bed-material 
size was greater than or equal to 60 millimeters (mm). Fur- 
thermore, the minimum value of K, that could be used in any 
instance was 0.70. 

The current (2011) HEC-18 equation with K, as recom- 
mended by the FHWA still has the form of equation 5, but the 
Mueller K, is computed by equation 11 for bed-material size 
conforming to current (2011) HEC-18 criteria (Richardson and 
Davis, 2001, p. 6.6). The criteria for applying K, require that 
the D,, particle size be greater than or equal to 2 mm and the 
D,, particle size be greater than or equal to 20 mm. The mini- 
mum value of K, under current (2011) criteria is 0.40. 


Recent Research of Pier-Scour Equations 


Pier-scour research has continued and new equations 
have been proposed that might lead to more reliable scour 
prediction than currently provided by conventional FHWA 
and HEC-18 equations. Recent (2003-05) pier-scour research 
in nonuniform sediment mixtures includes results reported 
by Molinas (2003) and Mueller and Wagner (2005). Molinas 
(2003) proposed two new methods that include equations 
similar in principle to those currently (2011) recommended 


by FHWA (equations 5 and 11). The first method proposed by 
Molinas (2003) includes the pier-scour equation: 


Vs 3 E 
p 064, 017 = K,K,K,K,Ky"” (13) 
n vs 
where 
K’, is aconstant equal to 1.21, in inch-pound 
units; 


y _ is the velocity-intensity ratio, dimensionless 


= V, = V on 
aa V ce = V im , un 


where 0 < y < 1 (this is the same ratio contained in equa- 
tion 8); and where all other terms are as previously defined. 
Also, velocity relations associated with equation 14 are: 


| 
Vice =K, Dery, (15) 
which is essentially the same as V”,, in equation 1, and 
V, =K,D,°Y", (16) 


which is the same as V,_ in equation 8 (where V, = 0.4V.,.) 
where 
K, is a constant equal to 12.0, in inch-pound 
units; 


K, is aconstant equal to 4.8, in inch-pound units. 


K, is then given by: 


D 
K,=1.25+3,|—* vIn (¥+0.5), 
Dsp 


for 0.0< K,< 1.0, and 0.0 < w< 1.0; and where all other terms 
are as previously defined. 


(17) 


The second method proposed by Molinas (2003) incor- 
porates K,computed by equation 17 into the HEC-18 equation 
(equation 5), with the addition of K, given by: 


V 0.45 
for V,>V, K,-(-%) (18) 


im? 


for VV < V. (18a) 


im? 


K=0.0, 


where K,is the correction factor for initiation of motion. 


K, from equation 17 and K, factor from equation 18 are incor- 
porated to produce the form of the HEC-18 equation proposed 
by Molinas (2003): 


2) 


0.65 
b 
y, = 2.0y, mx ck) ae (19) 


and where all other terms are as previously defined. The 
research by Molinas that resulted in equation 19 showed 
improvements among the laboratory data over the most 
revised version of the HEC-18 equation (equation 5 with K, by 
equation 6); however, that research was not supported by field 
data. 

Finally, Mueller and Wagner (2005) used field data from 
266 pier-scour measurements in the sand and coarser size 
range to show that the magnitude of overprediction by the 
HEC-18 equation without the inclusion of K, increased as the 
median size of the bed material increased. They proposed the 
development of a new methodology for determination of K, 
based on numerous combinations of variables analyzed. The 
best correlation resulted when K, was related to the pier width 
and median bed-material size according to: 


b, 0.19 
K, = 0.35 
Dayo 


where all terms are as previously defined. Mueller and Wagner 
(2005) determined that HEC-18 equation results using K, 
determined by equation 20 did not perform as well (statisti- 
cally) as the HEC-18 equation using K, determined by equa- 
tion 11 proposed by Mueller (1996). Still, the basis for the 
approach using equation 20 was considered justified because 
of previous research (Sheppard, 2004) that showed the ratio of 
pier width to sediment particle size to be an important scour 
mechanism. 


(20) 


Methods of Data Collection 


Data collected in the field for this study included infor- 
mation needed to measure pier scour along with bed material 
and ancillary data needed to compute scour using equations 
like those presented in the “Background” section of this report. 
Measured pier scour is a scour-depth estimate derived from 
various hydraulic measurements made at a bridge. Computed 
scour is a scour-depth prediction derived from variables mea- 
sured in the field and equations developed in previous studies. 
During 2001-07, pier-scour and related data were collected at 
bridges in Montana with a variety of stream sizes, streamflow 
conditions, pier configurations, and bed-material sizes and 
gradations. 

Limited-detail field procedures were used to investigate 
pier scour during this study. Measurements made at each site 
(tables 1, 2) included the following: 


1. Pier-scour depth derived from cross-section survey data. 
2. Bed-material characteristics. 


3. Approach velocity (V,) and approach depth of flow (y) at 
each pier. 


4. Water-surface elevation and slope at the bridge during 
high-flow conditions. 


5. Pier geometry and flow angle of attack on the pier (a). 


Introduction 9 


6. Background notes describing site-specific conditions. 


Upon returning from the field, equations | and 2 were 
used to compute the critical velocity to initiate scour of 
bed material at the pier (V’,.,), which was compared to the 
measured approach velocity (V,) to determine whether or not 
a valid real-time pier-scour measurement was obtained for a 
site. If a valid real-time pier-scour measurement was indicated, 
the information collected was used to determine measured 
scour and computed scour using equation 5. Where measured 
scour exceeded scour computed by equation 5, follow up data 
were collected to help explain the inconsistency. Follow up 
data-collection efforts included (1) measuring approach depth 
of flow at different magnitudes of streamflow to confirm the 
elevation of the reference bed surface (fig. 2), also called the 
concurrent ambient-bed surface (Landers and Mueller, 1996, 
p. 36-37; Mueller and Wagner, 2005, p. 8-9); (2) determining 
the potential effects of the thalweg on scour-depth estimation, 
(3) identifying the presence of substantial amounts of large, 
woody debris; (4) collecting additional bed-material samples; 
and (5) reassessing the flow angle of attack. 


Site Selection 


Bridge sites were selected to include a range of coarse 
bed-material characteristics, bridge-pier configurations, and 
stream characteristics such as bed-material size and gradation, 
pier width, and stream size that can affect pier scour. Several 
bed-material, pier, and site criteria were used to select suitable 
bridge sites for measurement of pier scour. 

Foremost was the bed-material criteria used to select 
sites. Only sites with noncohesive, nonuniform, coarse-bed 
material were considered, which is consistent with HEC-18 
criteria (Richardson and Davis, 2001), where the D., particle 
size is > 2 mm and the D,, particle size is > 20 mm. Most of 
the bridge sites investigated were located in stream reaches 
that would likely be classified as pool-riffle or plane-bed 
(Montgomery and Buffington, 1997). 

Criteria also were used to determine if piers associated 
with a site were suitable for scour measurement. When piers 
are aligned parallel with the flow, maximum scour typically 
occurs at the upstream nose of the pier (Mueller and Wagner, 
2005, p. 18). Bridge sites selected for measurement had piers 
with little or no flow angle of attack (a generally less than 
about 5 degrees). A few sites where a ranged from 5 to a maxi- 
mum of 10 degrees were selected because remnant scour holes 
with no substantial infilling indicated maximum scour depth 
was located near the upstream nose of the pier, and a few sites 
also were selected because flow became more aligned with the 
pier at higher stages of flow. Thus, scour measurements were 
made at the upstream end of the pier in all instances. 

Other criteria also were considered in the selection 
of sites. Sites selected had no unique hydraulic conditions 
that could produce pier scour by processes other than those 
accounted for in the pier-scour equations evaluated in this 
report. Furthermore, piers that had substantial accumulation 
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of large, woody debris or a tendency for debris accumulation 
based on historical information were not considered for scour 
measurements. Finally, sites were selected where contraction 
and abutment scour did not contribute to pier scour because 
(1) all piers were beyond the influence of abutments that may 
have projected into the main-channel flow, (2) the site had 

an incised channel with little or no flood plain or lacked any 
flow contraction at the bridge, (3) bridge abutments did not 
encroach within the bankfull width of the stream, or (4) runoff 
conditions at the site were insufficient to cause measurable 
contraction or abutment scour effects. 


Pier-Scour Depth 


Measured pier-scour depth (y",) is the difference between 
a measured reference surface of the streambed and the lowest 
measured location within a pier-scour hole. Scour measure- 
ments primarily consist of cross-section surveys, but also may 
include longitudinal profile surveys and continuous measure- 
ment of scour using fixed instrumentation. At all sites, scour 
measurements were conducted during intermittent site visits 
using portable instrumentation. At some sites, scour mea- 
surements also were made using instrumentation fixed at a 
particular location (Mueller and Landers, 2000 and Mueller 
and Wagner, 2005). 

Scour measurements obtained during high-flow condi- 
tions at each bridge fell into one of four categories: 


1. One measurement at a single pier. 
2. One measurement at multiple piers. 
3. Multiple measurements at a single pier. 


4. Multiple measurements at multiple piers. 


Most measurements were in category | or 2. 

The most common method to measure scour is to 
conduct cross-section surveys. For this study most cross- 
section surveys involved sounding with a cable-suspended 
Columbus-type (or C-type) weight (Jarrett and Boyle, 1986) 
deployed from the upstream side of the bridge deck. A few 
scour measurements were obtained by probing the streambed 
with a graduated metal rod when the bridge deck was close to 
the streambed surface and the streamflow depth was less than 
about 3 or 4 ft. 

In addition to these more-conventional methods, 
hydroacoustic devices also were used. A portable chart- 
recording echo sounder was used at several sites to survey 
cross sections. Two devices used with the echo sounder 
included a knee board with a transom-mounted transducer fit- 
ted through the board and a transducer mounted to the bottom 
of a Columbus weight (Mueller and Landers, 2000, p. 37-51). 
The echo sounder and transducer worked well in the stream 
channel, some distance away from the pier, but scour mea- 
surements at piers during high-flow conditions generally 
were problematic because (1) turbulence and drawdown near 
the pier made the knee board pitch, heave, or roll, so the 
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transducer was not submerged beneath the water surface or 
was out of vertical alignment at times; (2) air entrainment 
interfered with acoustic transmission; or (3) high surface 
velocities swept the kneeboard away from the pier. 

A boat equipped with a portable survey-grade echo 
sounder and global positioning system (GPS) was used in 
several of the large, deep rivers during moderate streamflow 
conditions to obtain channel-bottom data near the pier and 
longitudinally upstream and downstream from the bridge cross 
section. These data were used with other high-flow cross- 
section surveys to assess scour depth and to evaluate remnant 
scour. The prevalence of clear-water scour conditions result- 
ing in no substantial infilling allowed for measurements of 
remnant scour that would not have been valid under live-bed 
conditions. The boat and echo sounder were used at two large 
river sites exceeding 300 ft in flow width (sites 51 and 57, 
table 1), where the channel bottom was surveyed for several 
hundred feet upstream and downstream from the piers to aid in 
determining the reference bed surface longitudinally. The boat 
and echo sounder also were used at site 57 to map the pier- 
scour hole, where the flow angle of attack on three out of four 
piers was 10 degrees. At this site, soundings were conducted 
at depths in excess of 10 ft along the sides and along the 
upstream and downstream ends of the piers under moderate 
streamflow conditions to confirm that the deepest scour was 
located at the upstream end of the piers. 

Fixed instrumentation used to measure scour included 
acoustic transducers and scour rods installed at 10 sites. 
Transducers were used at three sites (sites 21, 33, and 41) and 
consisted of one or more 6-degree beamwidth, 235-kilohertz 
devices, fastened to the pier near the upstream end to measure 
scour depth on a near-continuous basis (15-minute time step). 
Scour rods were used at seven sites (sites 4, 5, 30, 37, 40, 43, 
and 45) and typically consisted of 3-ft- to 6-ft-long smooth 
steel rods 1-inch in diameter that were encased most of their 
lengths in a 2-inch-diameter steel pipe fixed to the pier near 
the upstream end to measure only maximum scour during 
spring runoff. Elevations of the top of the rod and casing were 
surveyed and tied to a local vertical datum before and after 
runoff and were used with other information to determine 
scour depth. The top of the scour rod and outer casing pro- 
jected above the water surface, or the top of casing was sealed 
to prevent intrusion of sand and finer suspended sediment dur- 
ing submergence. Streambed material supporting scour rods 
generally exceeded the diameter of the rod. Scour rods were 
observed in operation during very high approach velocities 
(velocity at site 45 was 11.1 ft/s, table 2) with no indication 
of erratic rod movement or vibration that might cause drilling 
action that would induce additional scour. 

Cross sections, piers, and other features were described 
using a USGS convention, in which features are referenced in 
a downstream view from left to right. Therefore, the left bank 
was considered to be on the left-hand side of the observer 
when looking downstream in the direction of streamflow. The 
left pier at a site, therefore, was designated pier number | and 
piers to the right were assigned pier numbers 2, 3, and so on. 
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Table 2. Pier-scour data collected for this study. 


[All pier-scour data are for sites with nonuniform, coarse-bed material (predominantly gravel and cobble). Pier number: 1, 2, 3 are piers left to right looking 
downstream. The term “single” is used when bridges had one pier. Pier type: S, single; G, group; Debris effects: I, insigificant; M, moderate. Pier-scour 

depth from field measurements. Abbreviations: ft, feet; ft/s, feet per second; G,, gradation coefficient, dimensionless; mm, millimeters; D,, surface bed-material 
particle size where i percent of particles are finer by weight, in millimeters] 


Site Measure- Pier Pier Pier Nominal Pier 
number Site name ment Date Time nose pier width, b length, L 
(fig.1) number number type shape (ft) " (ft) 
Missouri River Basin 

1 Big Hole River east of Wise River 1 06/04/02 1100 2 Ny) Sharp 2.00 15.6 

1 Big Hole River east of Wise River 2 06/04/02 1045 3 Ss) Sharp 2.00 15.6 

1 Big Hole River east of Wise River 3 05/24/06 1810 3 S Sharp 2.00 15.6 
2 Big Hole River west of Divide 4 06/04/02 1245 1 S Sharp 4.00 36 
2 Big Hole River west of Divide 5 06/04/02 1230 2 S Sharp 4.00 36 
2 Big Hole River west of Divide 6 05/25/06 1100 2 Ss) Sharp 4.00 36 
3 Beaverhead River north of Dillon 7 05/31/06 1130 2 Ny) Sharp 2.67 28 
4 Big Hole near Melrose 8 05/24/06 1040 1 S Round 7.00 24 
5 Boulder River near Basin 9 06/10/06 1845 1 S Round 1.69 84 
6 Wisconsin Creek near Sheridan! 10 06/29/05 1350 1 G_ Round 1.02 14 
7 Gallatin River northwest of West Yellowstone 11 06/28/06 1700 Single S Sharp 2.00 32 
8 Big Hole River southwest of Twin Bridges 12 05/31/06 1440 1 S Sharp 3.67 26 
9 Jefferson River north of Twin Bridges 13 06/13/06 1505 1 N) Sharp 3.00 46 
9 Jefferson River north of Twin Bridges 14 06/13/06 1520 2 S Sharp 3.00 46 
9 Jefferson River north of Twin Bridges 15 06/13/06 1535 3 Ny) Sharp 3.00 46 
10 Madison River south of Cameron 16 06/01/03 1245 1 S Sharp 2.80 33) 
10 Madison River south of Cameron 17 06/16/05 1730 1 S Sharp 2.80 33 
10 Madison River south of Cameron 18 06/16/05 1745 2 S Sharp 2.80 33 
11 Boulder River at I-90, near Cardwell 19 05/19/05 1445 1 Ny) Sharp 3.00 40 
11 Boulder River at I-90, near Cardwell 20 06/06/05 1220 il S Sharp 3.00 40 
12 Boulder River near Cardwell 21 06/06/05 1045 1 S Sharp 3.00 37 
3) South Boulder River near Cardwell 2D 06/29/05 1015 Single S Square 1.60 26 
13 South Boulder River near Cardwell 23 06/11/06 1650 Single S Square 1.60 26 
14 South Willow Creek near Harrison 24 06/11/06 1840 Single G Round 1.00 3i7/ 
13 Gallatin River west of Bozeman 25 05/30/06 1820 Single Ny) Sharp 3.00 41 
15 Gallatin River west of Bozeman 26 06/07/06 2025 Single S Sharp 3.00 41 
16 Jefferson River west of Three Forks 2) 06/13/06 2005 2 S Sharp 3.30 32 
7) Gallatin River near Logan 28 05/24/05 1715 Ds S Sharp 2.80 24 
18 Jefferson River north of Three Forks 29 06/28/05 1215 2 S Sharp 3:33 41 
19 Madison River near Three Forks 30 05/23/05 1315 Single S Sharp 3.00 32 
19 Madison River near Three Forks 31 06/15/05 1345 Single N) Sharp 3.00 32 
20 Missouri River near Townsend 32 06/06/05 1630 1 S Sharp 3.00 46 
20 Missouri River near Townsend 33 06/06/05 1700 2 S Sharp 3.00 46 
20 Missouri River near Townsend 34 06/12/06 1920 1 S Sharp 3.00 46 
20 Missouri River near Townsend 35 06/12/06 1935 2 S Sharp 3.00 46 
21 Sun River north of Augusta 36 06/12/02 1345 1 S Sharp 3.50 28 
21 Sun River north of Augusta 37 06/12/02 1415 2 S Sharp 3.50 28 
21 Sun River north of Augusta 38 05/25/06 1545 1 Ss) Sharp 3.50 28 
22 Marias River near Shelby 39 06/11/02 1015 1 N) Sharp 6.007 27 


Introduction 17 


Flow angle Approach Approach depth Debris Pier-scour Accuracy 
of attack,a velocity, V. of flow, y, effects le D,, D,, D,, 6, depth, y’, of y’, 
(degrees) (ft/s) (ft) (ft) (ft) 
Missouri River Basin 

0 7.50 751 I 40.5 102 176 263 2.09 1.14 0.30 
0 7.89 741 I 40.5 102 176 263 2.09 1.44 30 
0 7.79 7.09 I 40.5 102 176 263 2.09 1.12 30 
0 6.09 4.54 I 23M 71.3 155 190 2.56 1.36 30 
0 7.21 4.94 I 23.7 71.3 155 190 2.56 76 30 
0 6.65 4.76 I DBM) 71.3 155 190 2.56 1.35 30 
0 3.62 1.62 I 20.3 35.8 39:2 80.3 1.71 98 30 
0 5,115 4.10 I 28.2 59.5 104 157 1.92 D8) 20 
0 4.82 2.60 [ 29.1 58.0 111 163 1.96 32 10 
0 S27) 1.40 I 232) 47.8 78.8 104 1.84 40 20 
0 7.63 1.82 I 32.7 64.7 117 167 1.89 1.58 .20 
0 4.58 525 M 20.8 38.5 60.9 80.1 1.71 1.27 40 
10 5.98 Tadd I 28.8 60.7 114 148 1.99 1.35 40 
10 7.21 6.02 I 28.8 60.7 114 148 1,99) 233) 30 
8 6.80 7.83 [ 28.8 60.7 114 148 1.99 2.06 30 
0 503) 3.01 I DOG) 79.8 149 253 2223 .86 30 
0 5.64 3.22 I 29.9 79.8 149 253 2.23 78 30 
0 6.09 322 I DY) 79.8 149 53} DD a) 30 
0 4.00 4.73 I 2.58 17.1 44.1 82.9 4.14 1.37 30 
0 5.02 5.67 I 2.58 7 44.1 82.9 4.14 2.03 30 
0 4.45 4.71 Ne 2.58 17.1 44.1 82.9 4.14 71 30 
0 5.78 80 I DO 3) 35 111 155 1.95 40 20 
0 6.38 2.24 I 29.3 55.2 111 155 1.95 .64 30 
0 4.57 1.85 I 8.76 38.3 WS 109 Poi .60 20 
0 6.45 3.61 I 45.5 70.4 109 143 1.55 .68 30 
0 8.33 4.49 I 45.5 70.4 109 143 1.55 125) 30 
0 6.25 15:9 [ 5.91 22.3 57.0 89.6 3.11 2.08 40 
0 6.24 5.80 I 17.8 33} 57.4 76.0 1.80 1.70 30 
0 3.87 8.47 I 8.85 19.1 31.1 47.5 1.87 .90 30 
0 4.67 PATE) I 19.7 29.8 43.2 56.8 1.48 1.47 20 
0 5.22 1.76 I 19.7 29.8 43.2 56.8 1.48 1.10 .20 
0 5.83 8.47 I 11.9 29,7) 64.1 87.4 PSD, 1.23 50 
0 4.74 8.18 M 11.9 29.7 64.1 87.4 2.32 1.43 50 
0 7.84 10.3 I 11.9 DOT 64.1 87.4 230) DMO 50 
0 6.83 10.0 M 11.9 29.7 64.1 87.4 2.32 2.99 50 
0 5.74 3.60 I 26.8 57.4 111 156 2.03 1.50 20 
0 4.77 2.83 I 26.8 57.4 111 156 2.03 1.14 .20 
0 5.18 3,112 I 26.8 57.4 111 156 2.03 1.38 20 
0 9.56 14.7 I 7.36 22.7 44.0 59:3 2.44 2:97 50 
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Table 2. Pier-scour data collected for this study.—Continued 


[All pier-scour data are for sites with nonuniform, coarse-bed material (predominantly gravel and cobble). Pier number: 1, 2, 3 are piers left to right looking 
downstream. The term “single” is used when bridges had one pier. Pier type: S, single; G, group; Debris effects: I, insigificant; M, moderate. Pier-scour 

depth from field measurements. Abbreviations: ft, feet; ft/s, feet per second; G,, gradation coefficient, dimensionless; mm, millimeters; D,, surface bed-material 
particle size where i percent of particles are finer by weight, in millimeters] 


Site Measure- Pier Pier Pier Nominal Pier 
number Site name ment Date Time nose pier width, b length, L 
(fig.1) number number type shape (ft) " (ft) 
Missouri River Basin—Continued 
23 Smith River southeast of Ulm 40 06/08/05 0820 1 S Sharp 3.20 31 
23 Smith River southeast of Ulm 41 06/14/05 1115 1 S) Sharp 3.20 31 
24 Smith River at old truss bridge, southeast of 42 06/08/05 1415 2 S Round 3:93 16 
Ulm 
25 Belt Creek 7 miles south of Belt 43 06/13/05 1330 1 S Sharp 3.50 42.5 
26 Belt Creek 6 miles south of Belt 44 06/10/05 1445 1 S Sharp 3.50 42.5 
DT Belt Creek 2 miles southeast of Belt 45 06/13/05 1430 Single S Sharp 2.00 33.5) 
Yellowstone River Basin 
28 Yellowstone River at Emigrant 46 06/22/05 1630 2) Ss) Sharp 4.00 DS 
29 Shields River near Wilsall 47 06/30/05 1430 N) Square 2.50 25 
30 Mill Creek near Pray 48 06/22/05 1230 Single S Square 2.00 21 
30 Mill Creek near Pray 49 06/14/06 1600 Single S Square 2.00 21 
31 Yellowstone River near Pray 50 06/23/05 1000 | S Sharp 3.00 27 
31 Yellowstone River near Pray 51 06/23/05 0930 3 S Sharp 3.00 27 
31 Yellowstone River near Pray Sw 06/08/06 1400 1 S Sharp 3.00 2H] 
31 Yellowstone River near Pray 53 06/08/06 1420 3 Ny) Sharp 3.00 27 
ay Yellowstone River near Pine Creek 54 06/05/02 1145 1 S Round 3.50) 26 
32 Yellowstone River near Pine Creek oP) 06/05/02 1230 2 Ss Round 3.50 26 
33 Yellowstone River south of Livingston 56 06/10/06 1445 1 S Round 10.0 Ds 
33 Yellowstone River south of Livingston ay 06/14/06 1145 2 S Round 10.0 25 
34 Boulder River near McLeod 58 06/08/06 2015 1 S Sharp 3.00 22 
39 Boulder River south of Big Timber 59 06/08/06 1745 2 Ny) Sharp 3.60 34 
36 Boulder River at I-90 near Big Timber 60 06/09/06 1045 2 Ss) Sharp 3.60 44 
Columbia River Basin 
37) Little Blackfoot River near Avon 61 06/10/02 1500 1 G Round 1.69 46 
37 Little Blackfoot River near Avon 62 06/07/05 1135 1 G Round 1.69 46 
a7 Little Blackfoot River near Avon 63 06/08/07 1100 1 G Round 1.69 46 
38 Blackfoot River (old bridge) west of Lincoln 64 06/10/02. 1215 Single S Sharp 3.00 21 
38 Blackfoot River (old bridge) west of Lincoln 65 06/19/02 1445 Single S Sharp 3.00 21 
39 Blackfoot River (new bridge) west of Lincoln 66 05/13/05 1745 Single N) Sharp 2.00 42 
39 Blackfoot River (new bridge) west of Lincoln 67 06/05/06 1845 Single S Sharp 2.00 42 
39 Blackfoot River (new bridge) west of Lincoln 68 06/08/07 1745 Single S Sharp 2.00 42 
40 Little Blackfoot River near Garrison 69 06/11/06 1330 Single S Sharp 2.00 46 
40 Little Blackfoot River near Garrison 70 06/08/07 1320 Single S Sharp 2.00 46 
41 Clark Fork near Gold Creek 71 06/07/07 0745 2 S Round 5.00 31 
42 North Fork Blackfoot River west of Lincoln 72 05/30/02 1115 1 N) Sharp 3.00 42 
42 North Fork Blackfoot River west of Lincoln 73 05/30/02 1130 2 S Sharp 3.00 42 
42 North Fork Blackfoot River west of Lincoln 74 05/17/05 1400 2 S Sharp 3.00 42 
43 North Fork Blackfoot River near Ovando 75 05/20/05 1330 1 G Round 1.69 34 
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Flow angle Approach Approach depth Debris Pier-scour Accuracy 
of attack,a velocity, V. of flow, y, effects i D,, D,, D,, 6, depth, y’, of y’, 
(degrees) (ft/s) (ft) (ft) (ft) 
Missouri River Basin—Continued 

P 6.18 5.08 I 6.67 13.6 26.7 41.5 2.00 2.92 0.30 
5 5.89 4.57 I 6.67 13.6 26.7 41.5 2.00 2.24 30 
0 6.59 7.88 I 19.8 39.2 86.0 123 2.08 85 30 
57 6.88 M 15.1 44.6 88.8 125 2.42 2a) 30 

6.45 2.80 I 14.5 38.0 75.1 105 2.28 1.85 50 

0 6.11 2.86 I 23.1 39,5 76.3 109 1.82 1.14 20 

Yellowstone River Basin 
0 6.84 9.69 I 14.6 Bow 67.6 94.4 ANS 2.81 0.50 
0 5.04 1.37 I 24.8 50.6 88.0 115 1.88 AS .20 
0 6.10 DS2) I BORD, 67.4 130 164 2.01 wD) 20 
0 5.73 2.12 I 32.2 67.4 130 164 2.01 53 20 
0 6.48 6.32 I 41.5 OP. II 133 165 1.79 1.67 30 
0 7.61 8.97 M 41.5 92.1 133 165 1.79 1.40 40 
0 6.66 7.10 I 41.5 92.1 133 165 1.79 1.20 30 
0 8.13 9.77 M 41.5 92.1 133 165 1.79 1.63 40 
0 6.85 10.3 M 10.9 54.5 114 174 3) 98} Pao) 50 
0 8.25 11.5 I 10.9 54.5 114 174 3,23 3.05 40 
0 10.5 8.75 I 24.1 355 88.8 117 1.92 6.60 20 
0 5.44 6.00 I 24.1 53:5 88.8 117 1.92 3.13 40 
0 7.64 4.72 I 319) 55.6 108 128 1.84 80 A0 
0 ek? 8.71 I 43.2 76.6 125 171 1.71 99 30 
0 ee, 2.44 I 46.9 84.6 136 196 1.70 90 30 
Columbia River Basin 

0 3.62 2.48 I 13.2 DOD) 352 TID) 2.06 52 0.20 
0 6.00 1,72 I 13.2 29.2 55.2 TIO 2.06 1.08 .20 
0 5.06 3.18 I 132 DOD Sy TID 2.06 1.47 20 
0 3.03 7.80 I 1.68 9.63 26.8 38.2 4.00 1.00 30 
0 3,112 8.60 I 1.68 9.63 26.8 332 4.00 1.41 30 
0 3.04 5.24 I 1.68 9.63 26.8 38.2 4.00 .76 30 
0 3.78 4.76 I 1.68 9.63 26.8 38.2 4.00 84 30 
0 4.08 5.04 I 1.68 9.63 26.8 38.2 4.00 1.16 30 
0 4.87 2.41 I 14.6 57.4 98.1 a7) DS) .60 20 
0 5.98 2.70 I 14.6 57.4 98.1 127 2.59 75 .20 
0 4.22 5.58 I 8.62 30.9 79.8 159 3.04 2.09 20 
0 6.77 3.94 I 22.4 36.6 60.6 82.2 1.65 1.46 .20 
0 8.30 3.98 I 22.4 36.6 60.6 82.2 1.65 1.71 30 
0 5.81 2.84 I 22.4 36.6 60.6 82.2 1.65 56 30 
0 535 2.06 I MSY) 50.1 SB 128 1.90 AS 20 
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Table 2. Pier-scour data collected for this study.—Continued 


[All pier-scour data are for sites with nonuniform, coarse-bed material (predominantly gravel and cobble). Pier number: 1, 2, 3 are piers left to right looking 
downstream. The term “single” is used when bridges had one pier. Pier type: S, single; G, group; Debris effects: I, insigificant; M, moderate. Pier-scour 

depth from field measurements. Abbreviations: ft, feet; ft/s, feet per second; G,, gradation coefficient, dimensionless; mm, millimeters; D,, surface bed-material 
particle size where i percent of particles are finer by weight, in millimeters] 


Site Measure- Pier Pier Pier Nominal Pier 
number Site name ment Date Time nose pier width, b length, L 
(fig.1) number number type shape (ft) " (ft) 
Columbia River Basin—Continued 

44 Clark Fork south of Drummond 76 06/02/05 1210 2 S) Sharp 3.00 25 

44 Clark Fork south of Drummond 77 06/04/05 0945 2 Ny) Sharp 3.00 25 

45 Lost Horse Creek north of Darby 78 05/16/06 1015 1 G Round eS yi 46 

46 Clark Fork near Drummond 79 06/15/06 1045 2 S Round 1.67 26 

47 Clark Fork at Turah Bridge 80 06/05/01 1230 3 S Sharp 3.00 26 

48 Swan River east of Ferndale 81 05/27/06 1035 2 S Sharp 3.00 32 

48 Swan River east of Ferndale 82 05/27/06 1055 3 Ss) Sharp 3.00 3) 

49 Clark Fork at Missoula 83 05/23/02 0730 1 N) Sharp 3.25 68 

50 Clark Fork near Alberton 84 05/31/02 0830 2 S Sharp 3.00 26 

51 Clark Fork near Superior 85 06/01/06 1010 1 S Sharp 5.00 34 

51 Clark Fork near Superior 86 06/01/06 1055 2 S Sharp 5.00 34 

52 St. Regis River at Ward Creek 87 05/30/02 1930 2 S Square 1,35 25 

33 St. Regis River west of St. Regis 88 05/30/02 1745 Single S Sharp 3.00 29 

54 St. Regis River at St. Regis 89 05/30/02. 1515 Single N) Sharp 2.50 33 

55 Flathead River near Perma 90 06/22/06 1550 1 S Sharp 3.80 28 

55 Flathead River near Perma 91 06/22/06 1615 2 S Sharp 3.60 28 

55 Flathead River near Perma 92 06/22/06 1650 3 S Sharp 3.50 28 

56 Clark Fork above Flathead River, near 93 06/06/06 1215 2 S Sharp 3.50 34 
Paradise 

5 Clark Fork below Flathead River, near 94 06/02/06 1630 1 S Sharp 4.33 38 
Paradise 

57 Clark Fork below Flathead River, near 95 06/02/06 1615 2 S Sharp 4.33 38 
Paradise 

57 Clark Fork below Flathead River, near 96 06/02/06 1530 3 Ss) Sharp 4.33 38 
Paradise 

57 Clark Fork below Flathead River, near 97 06/02/06 1500 4 S Sharp 4.33 38 
Paradise 

37 Clark Fork below Flathead River, near 98 06/22/06 1245 2 S Sharp 4.33 38 
Paradise 

57 Clark Fork below Flathead River, near 99 06/22/06 1315 3 S Sharp 4.33 38 
Paradise 

58 Kootenai River near Libby 100 05/26/06 1430 es) S Sharp 4.00 34 

58 Kootenai River near Libby 101 06/21/06 1150 3 S Sharp 4.00 34 

59 Yaak River near Troy 102 05/26/06 1130 1 S Sharp 3.50 48 

59 Yaak River near Troy 103 05/26/06 1045 2 S Sharp 3.50 48 


'Site is abandoned railroad bridge on abandoned railroad grade. 


*Pier width based on exposed footing width of 6 ft. 
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Flow angle Approach Approach depth Debris Pier-scour Accuracy 
of attack,a velocity, V, of flow, y, effects ie D,, D,, D,, 6, depth, y’, of y’, 
(degrees) (ft/s) (ft) (ft) (ft) 
Columbia River Basin—Continued 
0 7.46 6.63 I 5.46 34.4 86.4 118 3.98 2.07 0.30 
0 5.79 ST I 5.46 34.4 86.4 118 3.98 1.53 30 
0 11.1 PAE) I 24.1 73.0 122 168 2.26 oS 20 
0 4.12 3.80 I 15.8 27.8 43.6 62.0 1.66 .70 30 
0 6.97 6.30 I 20.9 49.3 87.1 122 2.04 1.94 30 
0 7.09 7.36 I 34.8 129 115 146 1.82 1.15 30 
0 5.84 5.00 I 34.8 129 115 146 1.82 80 30 
0 8.14 15.4 I 28.5 67.3 123 168 2.08 195 40 
0 7.02 12.4 M 18.1 46.1 95.4 140 2.30 Dll 50 
0 5.98 9.76 I 20.0 53.4 105 144 2.29 2.10 50 
0 6.43 8.96 M 20.0 53.4 105 144 229 3.54 50 
0 9.37 7.30 [ 22.7 50.7 88.6 121 1.97 .90 30 
0 8.35 6.60 I D2) 41.0 69.6 86.4 1.60 1.58 30 
0 8.04 6.53 I PAT 24.5 50.0 84.8 1.99 137 40 
0 5.42 17.8 M 16.4 27.6 54.7 81.6 1.82 2.81 50 
0 4.90 19.4 M 16.4 27.6 54.7 81.6 1.82 1,92 50 
0 4.51 10.9 M 16.4 27.6 54.7 81.6 1.82 1.90 50 
0 6.34 10.7 I 16.2 54.9 115 159 2.67 2.29 30 
5 5.65 14.9 I 18.7 31.3 Sil 69.4 1.65 2.66 50 
10 6.73 13.5 M 18.7 31.3 51.1 69.4 1.65 7.00 50 
10 7.38 12.5 I 18.7 31.3 51.1 69.4 1.65 4.24 50 
10 6.03 17.6 M 18.7 31,3 51.1 69.4 1.65 6.43 50 
10 139) 14.8 M 18.7 31.3 51.1 69.4 1.65 6.50 50 
10 eke] 13.9 I 18.7 31.3 51.1 69.4 1.65 5.06 50 
0 6.91 6.93 I 27.4 58.0 93.0 126 1.84 1.87 30 
0 9.78 8.03 [ 27.4 58.0 93.0 126 1.84 2.47 50 
0 6.97 4.50 I 46.1 108 189 265 2.02 1.00 30 
0 6.51 3.95 I 46.1 108 189 265 2.02 1.05 50 
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The pier-numbering system used in the study differs from 
MDT pier numbers indicated on design drawings because 
bents or abutments at each end of the bridge structure included 
in the MDT pier-numbering system were excluded from the 
USGS pier-numbering system. The USGS pier-numbering 
system also was applied to sites operated and maintained by 
county and local entities that frequently had no engineering 
drawings with pier-numbering information. 

Cross-section data were used to determine the reference 
bed surface (fig. 2) from which the lowest measured eleva- 
tion of the scour hole was subtracted to determine measured 
or observed scour. The reference bed surface was sketched 
(Mueller and Wagner, 2005, p. 8) using cross-section data 
obtained during high-flow conditions and used to determine 
pier-scour depth (fig. 3). Surveys to establish the reference 
bed surface mostly were made along the upstream bridge face 
by using a manually operated four-wheel-base bridge crane 
and B-reel setup with a cable-suspended sounding weight 
deployed from the bridge deck (Rantz and others, 1982; 
Turnipseed and Sauer, 2010). Cross-section surveys varied in 
detail; complete cross sections were surveyed at a majority 
of sites. At a few sites, cross-section data were only obtained 
close to the pier, but far enough away (laterally) from the pier 
to accurately define the reference bed surface (figs. 2 and 3). 
In all instances, detailed soundings were conducted near the 
pier to define the scour hole. Measurements occasionally were 
obtained longitudinally along the streambed for some distance 
upstream and downstream along the major axis of the pier to 
help define the reference bed surface. 
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Surveys at most bridges were conducted on at least two 
separate dates—once during high-flow conditions (typically 
in May or June) to obtain pier-scour measurements, and once 
during low flow to assess channel-geometry changes, if any, 
attributed to general scour (Richardson and Davis, 2001, 

p. 5.1-5.18) or channel deposition relative to the high-flow 
measurement. Additionally, low-flow measurements were ana- 
lyzed in combination with high-flow measurements to detect 
the extent of scour-hole infilling, to indicate whether clear- 
water or live-bed scour conditions prevailed at a site, and to 
assess the relative channel stability at the bridge cross section 
over a range of flow conditions. 

Whether the total clear-water scour for a particular mea- 
surement resulted from a single scour event or multiple scour 
events was not addressed in this study. Cross-section measure- 
ments made during low-flow periods indicated remnant scour 
holes at most sites, so the total scour depth measured under 
high-flow conditions at a site could be the result of a single 
scour event or multiple scour events perhaps over consecu- 
tive below-average runoff years. As with similar studies using 
limited-detail field procedures to measure pier-scour depth, the 
assumption was made that scour depth produced by a single 
clear-water scour event could be the same as scour depth pro- 
duced by multiple clear-water events of similar flow magni- 
tude. Therefore, pier-scour depths reported in this study were 
attributed to high-flow conditions of similar magnitude leading 
up to and during the time scour measurements were made. 
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Figure 3. Surveyed cross section at upstream side of bridge showing pier-scour holes and reference bed surface 


used to determine pier-scour depth. 


Bed-Material Characteristics 


Bed-material sizes used in scour-prediction equations 
were determined from the surface bed material that generally 
reflected the armor layer. The potential interaction between 
the surface and shallow-subsurface bed-material layers (fig. 4) 
with respect to pier scour and armoring also was of interest 
for this study. Therefore, shallow-subsurface bed material was 
sampled and data were compared to surface bed-material data 
for each site. Surface and shallow-subsurface bed-material 
measurements were obtained by pebble counts and sieve 
analyses. Bed-material data were used to construct particle- 
size distribution curves, which were used to conduct various 
coarse-bed pier-scour analyses. 

Bed material was sampled at each site to characterize the 
average streambed material near the bridge, with an emphasis 
on the channel area leading into the bridge opening. Field 
observations indicated that coarse bed-material particle size 
and gradation varied at least as much vertically as areally near 
the bridge. Thus, surface bed-material samples were collected 
and analyzed to define bed-material characteristics in the same 
manner as traditional pier-scour studies. However, subsur- 
face bed material typically was sampled at relatively shallow 
depths less than | ft beneath the streambed surface. Both sam- 
ples generally were collected before or after the runoff season 
to more easily access and visually inspect the main-channel 
streambed. Samples were collected with square-hole templates 
and sieves to define the intermediate or b-axis (Bunte and Abt, 
2001, p. 20-27). 

The surface layer of the bed material was sampled by 
pebble counts (Wolman, 1954) where material was predomi- 
nantly fine gravel and coarser particles (4 mm and larger, 
table 3), and samples were sieved when the composition 
included sand and very fine gravel. Particles were measured in 
the field by pebble counts using the US SAH-97 gravelometer 
developed by the Federal Interagency Sedimentation Project 
(FISP) and described by Davis (2005). Pebble counts typically 
were limited to a single sample set of at least 100 particles 
obtained near the bridge opening. Although only a single sam- 
ple was obtained, considerable time was spent determining the 
optimal location to be sampled near the bridge. Samples were 
collected from the streambed in the main channel whenever 
possible. Locations were avoided where loose deposits of sedi- 
ment or the presence of finer material might indicate greater 
sediment mobility than was likely for the average bed-material 
composition. For flood-bar deposits, the mid-bar region was 
sampled thus avoiding the coarser and finer material typically 
deposited at the upstream and downstream ends of the bar 
(respectively). Multiple sets of pebble counts were conducted 
for about 10 percent of the bridge sites (table 4) used in the 
study. Some variation in particle size between the datasets at 
each site is indicated (standard deviation for D., table 4), but 
the relative magnitude of each standard deviation is small in 
most instances when compared to the corresponding mean 
particle size at each percentile (D,) shown. 
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Shallow-subsurface bed material was sampled within 
the zone typically considered part of the active layer or the 
upper layer of the streambed susceptible to sediment trans- 
port depending on hydraulic conditions (Molinas, 2000, p. 
33-37). Shallow-subsurface bed material was excavated by 
first removing the surface material or armor layer. Over- 
all depth of excavation from which the subsurface sample 
was collected was either equal to at least two times the D,, 
particle size, which generally describes the active-layer depth 
in coarse-bed streams (DeVries, 2002) or was about an 8- to 
10-inch minimum depth (whichever was greater). Bed material 
representative of the overall streambed was sampled whenever 
possible. Some samples were obtained along channel fringes 
or flood bars where the composition appeared similar to bed 
material located farther out in the stream, but where the com- 
bined effects of flow depth and velocity prevented sampling. 
A single sampling location typically was chosen, with other 
locations excavated and visually inspected to gain an indica- 
tion of average composition of the subsurface material. At a 
few sites, surface material was removed by merely scraping 
the surface material out of the way by hand. The majority of 
sites, however, required the use of a shovel to pry and exca- 
vate subsurface material. Many sites required the use of both a 
shovel and pry bar to obtain the sample. 

Surface material samples were used to characterize the 
subsurface material at four sites (site 31, site 37, site 38, and 
site 39; table 1) where excavation indicated that the size and 
gradation of the two layers were similar. Shallow-subsurface 
samples were sieved because sediment could not be character- 
ized by pebble counts alone. Sieved samples were air dried for 
several weeks or longer, hand sieved using standard 8-inch- 
diameter metal sieves for 11 screen sizes ranging from 2 to 
128 mm, and weighed. Samples were analyzed at the USGS 
Montana Water Science Center office in Helena, Mont (Guy, 
1969). Particle-size distribution curves were determined for 
the surface- and shallow-subsurface layers of each sample, and 
important bed-material size percentiles were interpolated from 
the curves. Percentiles included the D,,, D,,, Ds, D,, Dg; Dao, 
D,,, and D,, particle sizes, where the subscript i denotes the 
percentage of material finer (by weight) than the indicated par- 
ticle size (D). The D,, and D,. particle sizes were used to test 
the sensitivity of scour computations to vertical differences in 
bed-material size and to investigate how possible interactions 
between surface and shallow-subsurface streambed material 
are related to pier scour. 

Although the gradation and stratification of coarse 
streambed material can vary (Bunte and Abt, 2001, p. 129), 
all but three sites where pier-scour data were collected for this 
study had a higher fraction of coarse material on the surface 
than in the shallow subsurface. Such conditions are typical of 
armored streambeds, where enough coarse particles can accu- 
mulate to shield the entire bed surface (Lagasse and others, 
1991, p. 78-79). Although general channel stability was indi- 
cated by clear-water conditions and bed armoring upstream 
from the bridge at most sites, shallow-subsurface bed mate- 
rial was exposed at the piers when local scour took place. 
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Figure 4. Examples of differences in surface and shallow-subsurface bed material at sites 
where pier-scour measurements were collected for this study. A, Near-vertical perspective 

of the bed-material variation at site 15 (Gallatin River west of Bozeman). B, Plan view of bed- 
material variation at site 36 (Boulder River at I-90 near Big Timber). 


Table 3. Bed-material size classification used for this study 
(modified from Lane and others, 1947). 


[Symbols: =, greater than or equal to; <, less than] 


Particle-size range, Classification 


in millimeters 
>256 Boulders. 
128 to <256 Large cobbles. 
64 to <128 Small cobbles. 
32 to <64 Very coarse gravel. 
16 to <32 Coarse gravel. 
8 to <16 Medium gravel. 
4 to <8 Fine gravel. 
2 to <4 Very fine gravel. 
0.062 to <2 Very fine to very coarse sand. 
<0.062 Silt and clay. 


The shallow subsurface layer may, therefore, be important in 
determining appropriate bed-material variables to be used in 
the scour equations and was investigated accordingly. 

Based on the D,, particle size at each site, a Shields 
parameter (9) equal to 0.047 (Meyer-Peter and Miiller, 1948; 
Gessler, 1971) was used with equations | and 2. This value 
is consistent with a recent study (Mueller and others, 2005), 
in which a median value of @ equal to 0.040 was reported for 
45 gravel-bed streams and rivers. 


Ancillary Data 


Other data also were collected that were necessary to 
determine measured pier-scour depth and to compute scour. 
These data included approach velocity (V,) and approach 
depth of flow (y,) for all sites, water-surface elevation and 
slope, pier geometry, and flow angle of attack on the pier (a). 
Most of the ancillary data were collected during site visits, but 
some data at sites with fixed instrumentation were measured 
continuously. 

Approach velocity ideally represents the velocity just 
upstream from a pier, beyond the extent of the pier-scour 
hole and flow-acceleration zone (fig. 2). In the field, approach 
velocity typically is measured from a bridge deck during high- 
flow conditions. A vertical-axis mechanical current meter and 
USGS standard streamgaging procedures (Rantz and others, 
1982 ; Turnipseed and Sauer, 2010) were used to measure the 
velocity on either side of the piers. Velocity measurements 
adjacent to piers but outside the scour hole were made to 
approximate the velocity just upstream from the pier where 
velocity was difficult to measure. Approach velocity measure- 
ments made using a current meter were checked at one site 
(site 55, table 1) using a portable acoustic Doppler current 
profiler (ADCP). The specific ways in which approach veloc- 
ity was determined varied with the method used to measure 
pier scour at a site. 
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At sites where pier scour was measured using cross- 
section surveys during high-flow conditions, approach 
velocities were measured during the same site visit. Velocity 
was measured at several locations on either side of a pier 
to determine the approach velocity reported for the study. 
Multiple measurements of velocity were made to describe the 
variability of the approach velocities and any anomalies at a 
particular location. 

At sites where pier scour was measured with scour rods, 
approach velocities were determined using stage-velocity 
rating curves. At these sites, approach velocity and associated 
water stage were measured from the bridge deck during peri- 
odic site visits that included a range of moderate- to high-flow 
conditions. This limited number of measurements was used to 
develop a stage-velocity rating curve. The maximum seasonal 
water stage surveyed at the bridge, based on high-water marks, 
was then used in the stage-velocity rating curve to estimate 
maximum approach velocity by rating-curve extension. This 
method is consistent with findings that indicate maximum 
scour depth is closely correlated with maximum approach 
velocity (Miieller, 1996, p. 91). 

At sites where pier scour was measured with acoustic 
transducers mounted to piers, approach velocities were deter- 
mined based on periodic measurements of velocity, continuous 
measurement of stage, and development of stage-velocity rat- 
ing curves. The recorded water stage at the bridge correspond- 
ing to maximum scour determined from transducer soundings 
was used with the stage-velocity rating curve to estimate, by 
rating-curve extension, approach velocity associated with 
maximum pier scour. 

Approach depth of flow ideally represents the depth to 
stream bottom just upstream from the pier beyond the extent 
of the pier-scour hole and flow-acceleration zone (fig. 2). In 
the field, approach depth of flow typically is measured by 
sounding from a bridge deck during high-flow conditions. 
Bottom depth was sounded at several locations on either side 
of a pier to determine the approach depth of flow reported 
for the study. Multiple measurements of depth were made to 
describe the variability of the approach depth of flow and any 
anomalies at a particular location. Bottom depth was refer- 
enced to a water-surface elevation to determine the approach 
depth of flow. For high-flow scour measurements using cross- 
section surveys, the referenced water-surface elevation was the 
stage noted on the day and time when the bottom depth was 
sounded. For sites where pier scour was measured with scour 
rods, the referenced water-surface elevation was the maximum 
stage based on high-water marks surveyed after runoff. For 
sites where pier scour was measured with acoustic transduc- 
ers, the referenced water-surface elevation was the recorded 
stage at the time that maximum scour was indicated by the 
transducer. 

Water-surface elevation and slope during soundings 
were surveyed by trigonometric leveling with a total station 
and prism (Ghilani and Wolf, 2008) or by measuring from a 
reference mark down to the water surface. Vertical control for 
elevations and water-surface elevation measurements were 
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Table 4. Variation in surface bed-material particle size for multiple pebble counts made at bridge sites where pier-scour 


measurements were made for this study. 


[Abbreviations: D,, surface bed-material particle size where i percent of particles are finer by weight; mm, millimeters; n=, number of pebble counts in the 


dataset. Symbols: O,5 gradation coefficient] 


Site 1 Site 15 Site 17 
n=3 n=3 n=3 
Particle Mean value Standard Particle Mean value Standard Particle Mean value Standard 
‘ for D, deviation for D., : for D, deviation for D, : for D, deviation for D, 
size ene: ! size Soll : i size ced : : 
inmm inmm inmm inmm inmm inmm 
7 41 3.4 Di 46 4.8 Dig 18 0.7 
Dy, 73 6.8 Dy 58 3.0 Dig 27 1.3 
Di 102 7.3 D.. 70 5.8 Dy 33 22 
De, 176 15 Daa 109 19 De 57 2.4 
D, 179 14 Des 112 19 Ds 58 2.5 
Dy 209 8.4 Dy 127 22 Dyy 64 3.9 
Dag 263 4.9 Di 143 24 Dis 76 12 
Dy 309 29 Da, 179 13 Da, 102 19 
Mean o,= 2.09 Mean o,= 1.59 Mean o,= 1.80 


Standard deviation = 0.01 


Standard deviation = 0.09 


Standard deviation = 0.01 


Site 28 Site 36 Site 51 
n=4 n=3 n=3 
Particle Mean value Standard Particle Mean value Standard Particle Mean value Standard 
. for D, deviation for D, : for D, deviation for D, : for D, deviation for D, 
size a : i size eee i size ee : i 
inmm inmm inmm inmm inmm inmm 
- 15 3.0 D,, 47 18 D,, 20 4.0 
Dy 25 2.9 Ds 68 17 Ds 37 5.8 
Ds, 33 47 De 85 15 De 53 Sl 
Dea 68 21 oe 136 8.7 Di, 105 4.1 
D.. 69 22 D., 140 7 D., 108 3.3 
Dy 79 26 Dy 164 6.9 Dy 121 9 
DD. 94 32 Dag 196 26 Dag 144 9.1 
oe 110 42 ae 266 31 a 174 4.2 
Mean o,= 2.16 Mean o,=1.75 Mean o,= 2.31 


Standard deviation = 0.40 


Standard deviation = 0.25 


Standard deviation = 0.18 


referenced to a local arbitrary datum at some sites, whereas 
bench marks and other USGS gage reference marks were used 
where available. Water-surface slope is not an explanatory 
variable typically used in scour-prediction equations, but is 
needed to compute the shear stress on the streambed and to 
establish initial boundary conditions for site-specific hydraulic 
modeling studies. Therefore, water-surface slope was surveyed 
when the high-flow scour measurements were made to enable 
future research. 

Pier geometry was measured or obtained from engineer- 
ing drawings, and flow angle of attack on the pier (a) was 
determined visually in the field. The shape of the pier nose 
was recorded. At sites with tapered piers, nominal pier width 


(5) was computed as the average width in contact with flow at 
each surveyed water-surface elevation or measured stage. The 
flow angle of attack on the pier was determined during high- 
flow measurements. At some sites, the flow angle of attack 
also was assessed during low-flow conditions by observing the 
relative symmetry of the scour hole and the relative alignment 
of sediment deposited alongside and downstream from a pier. 
Relative pier scour (y’,/b) and relative bed-material size 
(b/D,,) were calculated for each measurement and used to 
relate measured pier-scour depth (y",), effective pier width (5), 
and median bed-material size (D,,). These variables normalize 
the data and provide dimensionless plots that indicate a rela- 
tion with pier-scour depth. Relative pier scour for study data 


was computed by using the effective pier width (b) according 
to: 
b=b,cosa+Lsina, (21) 
where 
L is the length of the pier measured 
perpendicular to the nominal pier width 
(L> 5), in feet; 
a. is the flow angle of attack on the pier, in 
degrees; 


and where all other terms are as previously defined. The effec- 
tive pier-width adjustment (where b > 6 , when a > 0 degrees) 
accounts for increased scour when flow is not parallel to the 
major axis of the pier. 

Collection of limited-detail measurements included 
recording background field notes and documenting hydrau- 
lic and scour conditions at each site. Typical observations 
included weather conditions, streamflow conditions, extent 
of any woody debris accumulated on piers, description of 
remnant scour holes, and bed-material characteristics and bed 
form in the stream reach near the bridge. For example, the 
bed-form correction factor (K,) was assigned a value of 1.1 
for all measurements, based in part on the coarse-bed material 
and plane-bed conditions recorded in the field. The quality 
rating of velocity measurements and soundings also was often 
described, for later use in assessing the suitability of measure- 
ments for determination of approach velocity and approach 
depth of flow. 

Water temperature can affect sediment transport and was 
therefore determined. Water temperature was either measured 
during a scour measurement with a calibrated thermometer, 
or estimated by analyzing real-time water temperature records 
reported for USGS streamflow-gaging stations near the mea- 
surement site. Based on an analysis of site measurements and 
water temperature records, water temperature was estimated to 
average about 50°F over the runoff periods for all sites. 


Pier Scour in Coarse-Bed Streams in 
Montana 


A total of 103 pier-scour measurements were made at 
59 bridges on coarse-bed streams in Montana during 2001-07 
for this study. Runoff conditions during the study mostly were 
limited to bankfull discharge approximating the 1.5- to 2-year 
recurrence-interval flood, which have a 67- and 50-percent 
chance, respectively, of happening in any year. Data collected 
for this study were compared to BSDMS scour measurements 
for coarse-bed streams and envelope curves based on previ- 
ous studies (Ettema, 1980; Melville, 1984; Chiew, 1984; and 
Mueller and Wagner, 2005). Data from this study also were 
compared with selected historical pier-scour measurements 
from Montana associated with peak streamflows of greater 
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recurrence interval (up to about 100 years) than pier-scour 
measurements collected for this study. Surface and shallow- 
subsurface bed-material data were collected for each site 
where pier-scour measurements were made to characterize the 
particle size and gradation. Bed-material characteristics in the 
two layers were used to compute scour by using the current 
(2011) version of the HEC-18 equation (equation 5), and the 
computed scour was compared to the measured scour. Rela- 
tions between surface and shallow-subsurface bed-material 
characteristics and pier scour also were analyzed. 


Hydrologic and Hydraulic Conditions 


Seasonal runoff in the mountain and foothill regions of 
the western one-half of Montana (fig. 1) typically results from 
spring snowmelt supplemented with May and June rainfall. 
In that region, summary data (McCarthy, 2005) indicate that 
runoff can begin in April, typically peaks in May or June, and 
recedes through the summer months to base-flow conditions 
in the fall. Isolated thunderstorms and broader-based general 
storms can cause occasional periods of high flow with vary- 
ing durations in the summer and fall. Streamflow conditions 
during this study generally were below average for the three 
major river basins in the study. Maximum daily mean flow 
during the spring runoff seasons was generally below aver- 
age, and spring-runoff volumes were below to near average 
(Attp://waterdata.usgs. gov/mt/nwis). Annual flood-frequency 
data (Parrett and Johnson, 2004) indicate that instantaneous 
peak discharges at most gaging stations near bridge sites 
where scour was measured approximated the 1.5- to 2-year 
recurrence interval during 2001-07 runoff. Lawlor (2004) 
found that the recurrence interval for bankfull discharge at 
41 streams in western Montana generally ranged from 1.0 to 
4.4 years with a median value of 1.5 years. Stream velocities 
and seasonal high-water marks measured at some ungaged 
sites in smaller drainage basins indicated that more intense 
runoff conditions, and hence, higher recurrence-interval peak 
streamflows occurred in these smaller basins compared to peak 
flows in larger drainage basins. Despite general below-average 
runoff conditions during this study, streamflow at most sites 
probably reached bankfull to near-bankfull conditions, which 
can produce notable sediment transport (Andrews, 1983) and 
scour. For example, the rate of increase in bedload discharge, 
based on 1999-2001 data at one study site (site 32, fig. 1 
and table 1), increased as streamflow approached bankfull 
discharge and then bedload discharge diminished beyond the 
bankfull stage (Holnbeck, 2005, p. 12-16). Real-time pier 
scour recorded at a downstream site (site 33, fig. 1 and table 1) 
by a fixed acoustic transducer between May—June 2006 was 
closely linked to changes in streamflow (fig. 5) for flows up 
to bankfull. The response of scour-hole depth to changes in 
streamflow at this site is complicated, however, by the sedi- 
ment supply from upstream that also increases during higher 
streamflows. For example, the scour hole is shown to be 
infilling (fig. 5) at times during the rising limb in early June, 
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Figure 5. Pier-scour depth and 
streamflow measured at site 33, 
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most likely as a result of bedload transport from upstream. The 
plateau that persisted for 3 days in early June (fig. 5) likely is 
due to woody debris deposited in the scour hole and later dis- 
lodged or to temporary armoring when small cobbles that can 
be transported (Holnbeck, 2005, p. 8-20) partially infilled the 
scour hole. Transducers at three other sites recorded similar 
responses of scour to changes in streamflow during the study 
period. 

The persistence of below- to near-average streamflow 
conditions during the study period likely produced repeated 
clear-water scour conditions at most study sites over succes- 
sive years. Therefore, infilling of scour holes from one runoff 
season to the next was probably minimal or only partial, 
allowing for the continued scouring of holes during succes- 
sive clear-water periods. Overall, runoff conditions during the 
study probably limited the opportunities to measure maxi- 
mum scour, but sites where measurements were successfully 
made reflect real-time scour conditions at piers in coarse-bed 
streams. 


Scour Data Collected 


Pier-scour measurements were made on various dates at 
59 bridge sites on coarse-bed streams. Of those measurements, 
103 satisfied criteria requiring the measured approach velocity 
(V,, table 2) exceed the critical velocity (V"_..,) to initiate scour 
of bed material at the piers. Incipient motion computations 
based on the D,, particle size indicated that 96 out of 103 pier- 
scour measurements were made under clear-water scour condi- 
tions. Two of the measurements were made under hydraulic 
conditions of supercritical flow, one measurement was made 
at or very near critical flow, and 100 measurements were made 
under subcritical flow. Comparisons of cross-section data 
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2-year recurrence interval peak streamflow 
approximately bankfull discharge 


Runoff 
hydrograph 
(15-minute 
time step) 


Maximum measured 
pier-scour 
depth =6.60feet Plateau likely due to 
debris or temporary 


scour-hole armoring 


May 2006 June 2006 


obtained during high-flow measurements with pre- or post- 
runoff cross-section data generally indicated a lack of substan- 
tial scour-hole infilling, which also supports the conclusion 
that clear-water conditions prevailed. Of the measurements 
made, 50 percent had an approach velocity (V) that equaled or 
exceeded 70 percent of the critical velocity (V,.,) for incipient 
motion of bed material upstream from the bridge, which might 
indicate that scour was measured very near the threshold 
between clear-water and live-bed scour (V'/V_., = 1.00) where 
maximum pier scour was shown in laboratory studies (Chiew, 
1984; Melville, 1984). 

The pier and stream characteristics associated with the 
103 measurements made for this study were variable but some 
general observations can be made about the measurements 
based on the individual measurements and summary statistics 
(tables 2 and 5). Despite reasonably high approach veloci- 
ties, measured pier scour was less than a couple of feet in 
depth for 70 percent of the measurements made at the coarse- 
bed sites. Particle size of sampled surface bed material was 
predominantly very coarse gravel based on the D., particle 
size. Surface bed material at one-half of the sites lacked 
substantial gradation, as indicated by gradation coefficients 
less than about 2.0. Drainage areas for most sites ranged from 
a few hundred square miles to a few thousand square miles 
in drainage area. About 60 percent of the sites were located 
on high-gradient streams based on water-surface slopes that 
exceeded 0.002 (Jarrett, 1984). The approach depth of flow 
typically was less than twice the effective pier width. This 
condition is probably common in high-gradient streams in 
the northern Rocky Mountains where flow depths tend to be 
shallower compared to flatter-gradient streams having similar 
streamflow magnitudes and bridge pier widths. There were no 
instances where pier-scour measurements involving relatively 


shallow approach flow depths in relation to effective pier 
width required later adjustment of computed pier-scour depth 
for wide-pier effect (Richardson and Davis, 2001, p. 6.7). All 
sites were subjected to hydraulic conditions where real-time 
scour was in progress at the time of measurement, as indicated 
by relative velocities (V/V...) in excess of 1.00 (table 5). 
Maximum relative pier scour from this study is less than 
values associated with previously published envelope curves, 
and therefore, data from this study do not redefine the upper 
magnitude of those envelope-curve relations. 

The effects of high-flow duration on pier scour were not 
investigated because most scour data for a site were collected 
in a single day. Continuous scour measurements made at three 
sites using fixed acoustic transducers, however, indicated that 
scour typically increased as streamflow increased. Pier scour 
at site 33, with a 10-ft wide pier (fig. 5, table 2), included 
remnant scour of about 5 ft, with scour depth fluctuating from 
a minimum of about 4 ft (some infilling) to a maximum of 
about 6.6 ft. Maximum pier scour at this site was attained for 
hydraulic conditions near the threshold between clear-water 
and live-bed scour (V/V, = 1.05). Transducer soundings indi- 
cate that a similar magnitude of maximum pier scour (about 
6.6 ft) was attained several times during May—June 2006 
(fig. 5) when flows ranged from less than to slightly greater 
than bankfull discharge. These continuous data from the 
transducers provide some insight on the site-specific effects of 
high-flow duration on pier scour. 

The effective pier-width adjustment for flow angle of 
attack (equation 21) was made for 11 of the 103 scour mea- 
surements. The adjustment was applied when a ranged from 
5 degrees to a maximum of 10 degrees (table 2) and when the 
flow angle of attack probably contributed to a small to moder- 
ate increase in scour depth at the upstream nose of the pier. 
When reporting relations between pier and streambed char- 
acteristics and scour depths, effective pier width was always 
used. 


Comparisons with Other Pier-Scour Data 


Pier-scour data collected for this study were compared 
to pier-scour data from other USGS studies that were used to 
develop pier-scour equations for coarse-bed streams. Histori- 
cal data were retrieved from the BSDMS described by Landers 
and others (1996). The 508 measurements in the BSDMS data- 
base collected between 1965 and 1998 by USGS were initially 
screened using the following criteria: 


1. Bed material was noncohesive and composed of 
sand and coarser particles. 


2. Data for particle-size distribution of the bed material 
were reported in the BSDMS database. 


3. Flow angle of attack on the pier (a) was reported 
as zero degrees. Sites with flow angles of attack 
exceeding zero degrees were not considered because 
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it was not known what adjustments might have been 
made to account for a effects on scour depth. 


4. Large, woody debris accumulation on piers was not 
indicated to be “substantial.” 


5. Measured approach velocity (V,) exceeded the criti- 
cal velocity to initiate scour of bed material at the 
pier (V"..,), where V/V"... > 1.00. 


6. All variables needed for computing scour depth were 
reported in the BSDMS database for the measure- 
ment or could be derived from other variables 
reported in BSDMS. 


The initial screening resulted in 272 pier-scour mea- 
surements that met the above criteria for highway bridges in 
15 States, referred to hereinafter as the BSDMS AI dataset. 

The measurements used for this study and measurements 
from the BSDMS All dataset were classified based on the 
size of bed material relative to pier width using the system 
described by Ettema (1980). In that classification, sediment 
size is considered: fine for b/D,, = 130; intermediate for 
130 > b/D,, 2 30; and coarse for 30 > b/D,, = 8. For b/D,, <8, 
sediment size is considered large and sediment typically is not 
eroded. 

A substantial portion (>75 percent) of the measure- 
ments in the BSDMS All dataset was made at sites with fine 
bed material based on Ettema’s classification system (fig. 6, 
table 6). In contrast, about two-thirds of pier-scour mea- 
surements made for this study were at sites with coarse bed 
material, and a small percentage (about 7 percent) of those 
measurements were at sites on the threshold between coarse 
and large bed material. 

All the data collected for this study and the BSDMS 
All dataset plot below the line reported by Raudkivi and 
Sutherland (1981) that showed the relation between relative 
pier-scour depth and relative bed-material size (fig. 6) near the 
clear-water and live-bed threshold (V/V. = 0.9). Data plot- 
ting below the line indicate that scour is not underpredicted. 

After the initial screening, the 272 measurements that 
make up the BSDMS All dataset were screened again for mea- 
surements where the D., and D,. particle sizes were equal to 
or greater than 2 mm and 20 mm (respectively) so K, could be 
less than 1.0. This screening reduced the BSDMS AI dataset 
from 272 to 90 measurements. All subsequent comparisons of 
BSDMS data either use the BSDMS All dataset for noncohe- 
sive bed material (n=272) or the final screened dataset (n=90), 
referred to hereinafter as the BSDMS Coarse dataset. The 
BSDMS Coarse dataset is comparable with data collected for 
this study based on FHWA criteria for D,, and D,. particle size 
that defines coarse bed material. 

In addition to using Ettema’s classification to compare 
relative bed material size (fig. 6), surface bed-material data 
collected for this study were compared with the BsDMS 
All and Coarse datasets (table 7) using the D., particle size 
and the bed-material size classification (table 3) to gener- 
ally characterize bed material size. Using very coarse gravel 
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EXPLANATION 


© Data collected for this study (7=103) 
x BSDMS All dataset (n=272) 


Variables 
n is the number of measurements 
be is pier-scour depth from field measurements, in feet 
b is effective width of the pier measured perpendicular to 


the approach flow, in feet 


D is the bed-material particle size where 50 percent of 
particles are finer by weight, in millimeters or feet 


RELATIVE PIER-SCOUR DEPTH (y’, /b), DIMENSIONLESS 


1 10 100 1,000 


RELATIVE BED-MATERIAL SIZE (4/Dgq), DIMENSIONLESS 


Figure 6. 
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Comparison of pier-scour measurements collected for this study and from the Bridge Scour Data Management System 


(BSDMS) All dataset. Relative bed-material size criteria from Ettema (1980). 


Table 6. Comparison of the relative bed-material size associated with pier-scour measurements made in this 
study with the Bridge Scour Data Management System (BSDMS) All dataset. 


[Symbols: >, greater than; >, greater than or equal to; <, less than; (b/D.,), relative bed-material size; n=, number of data values in 


dataset] 


Interval for relative bed-material Classification for 


Percent of data in size classification 


size (b/D,,), dimensionless bed-material size’ This study (n=103) BSDMS All dataset (n=272) 
b/D,, 2 130 Fine 0.00 76.5 
130 > b/D,, = 30 Intermediate 33.7 16.9 
30 > b/D,, 28 Coarse 59.2 6.6 
b/D,, <8 Large 7.1 .00 


'Bed-material size critera defined by Ettema (1980). 


as the size reference, about 67 percent of the measurements 
from this study are classified very coarse gravel or larger, only 
about 8 percent of the BSDMS All dataset is classified very 
coarse gravel or larger, and about 26 percent of the BSDMS 
Coarse dataset is classified very coarse gravel or larger. The 
BSDMS datasets do not include data for shallow-subsurface 
bed material, however, 54 percent of the shallow-subsurface 
measurements from this study are classified very coarse gravel 
or larger. None of the data from this study or the comparable 
BSDMS Coarse dataset (table 7) are ripple-forming sediments 
(D,, < 0.6 mm), whereas a substantial portion (about 32 per- 
cent) of the BSDMS All dataset are ripple-forming sediments. 


Ripple-forming sediments are more applicable to sand-bed 
streams, where fine noncohesive particle sizes and bed forms 
such as dunes are involved in the pier-scour process. 

Approach velocity (V.) is among the most important 
hydraulic variables that can affect pier-scour depth. The 
median approach velocity for measurements collected for this 
study was slightly lower than the median approach velocity for 
measurements in the BSDMS Coarse dataset but was higher 
than the measurements in the BSDMS AIl dataset (fig. 7). The 
approach velocity for the measurements in the BSDMS All 
dataset is more variable. 
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Table 7. Comparison of bed-material size classification associated with pier-scour measurements made for this study and 
measurements from the Bridge Scour Data Measurement System (BSDMS) Coarse and All datasets. 


[Abbreviations: mm, millimeter. Symbols: =, greater than or equal to; <, less than; D 
values] 


o median particle size of the sample; --, not applicable; n=, number of data 


Number of pier-scour measurements in each bed-material size classification’ 


Sod nintatialaira wlassilicatan Surface bed material Shallow-subsurface bed material’ 
This study BSDMS All dataset BSDMS Coarse dataset This study 
(n=103) (n=272) (n=90) (n=103) 

Small cobble 26 16 16 21 
Very coarse gravel 43 7 7 35 
Coarse gravel Di 35) 35) 3) 
Medium gravel 7 13 13 12 
Fine gravel 0 19 19 0 
Very fine gravel 0 0 

Sand (very fine to very coarse) 0 182 0 
Ripple-forming sediments? 0 88 -- 
Nonripple-forming sediments* 103 184 90 -- 


‘Based on D,, particle size of sample. 
°No data reported in BSDMS for shallow-subsurface bed material. 
°D,, < 0.6 mm. 


*D,, = 0.6 mm. 
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Figure 7. Comparison of approach velocities for pier-scour measurements collected for this study with the Bridge Scour Data 
Management System (BSDMS) Coarse and All datasets. 


Relative scour velocities (V,/V"..,) determined from 
measurements collected for this study also were compared 
with the BSDMS Coarse dataset (fig. 8). The threshold for 
initiating scour of bed material at the pier (V,/V".,. > 1.00) was 
exceeded to a greater degree by measurements in the BSDMS 
Coarse dataset. These higher relative scour values are likely 
due to the smaller bed-material size of measurements in the 
BSDMS Coarse dataset, compared to measurements collected 
for this study. Lower approach velocities can initiate scour at 
a pier when bed material is smaller. The median D,, and D,, 
bed-material size for data collected for this study are about 1.8 
and 1.6 times (respectively) larger than the D,, and D,, of the 
BSDMS Coarse dataset (fig. 9). 

In addition to comparing data used for this study to the 
BSDMS All and BSDMS Coarse datasets, data collected for 
this study also were compared to historical pier-scour mea- 
surements made in Montana by USGS during large floods. 
The measurements were not part of the BSDMS database. 
Pier-scour measurements made at sites 60, 61, and 62 (fig. 1, 
table 8) in 1996-97 were selected for this comparison. Nine 
measurements made at these three sites were for peak stream- 
flows approximating the 100-, 10-, and 5-year recurrence 
intervals, which are much greater than any recurrence intervals 
associated with measurements made for this study and provide 
an independent check on the envelope curves bounding the 
upper end of both the data collected for this study and BSDMS 
data combined. 
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While conducting research leading to the development of 
K,, Mueller (1996) introduced the concept of an idealized K;: 


' 


Idealized K, = —, (22) 
Vs 
where 

y', 1s pier-scour depth from field measurements, 
in feet; 

y, is pier scour computed by using the basic 
form of the HEC-18 equation without K, 
applied (equation 4), in feet; and 

K is dimensionless. 


4 


The idealized K, is the relative error between measured scour 
and computed scour. Idealized K, values were computed for 
pier-scour measurements collected for this study along with 
Mueller’s velocity-intensity ratios (equation 10). The resulting 
idealized K, and Mueller’s velocity-intensity ratios were then 
plotted in relation to each other for the pier-scour measure- 
ments collected for this study, the BSDMS Coarse dataset, and 
historical pier-scour measurements from sites 60-62 (fig. 10). 
The dashed line (fig. 10) defines an idealized K, according to 
HEC-18 criteria currently (2011) used (Richardson and Davis, 
2001), and represents the envelope curve under which all ide- 
alized K, values derived from field data should plot. 

Data collected for this study included lower Mueller 
velocity-intensity ratios than the BSDMS Coarse dataset and 
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Figure 8. Comparison of relative scour velocities for pier-scour measurements collected for this study and from the Bridge Scour Data 


Management System (BSDMS) Coarse dataset. 
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Figure 9. Comparison of bed-material particle size (D,, and D,,) for pier-scour measurements collected for this study with the Bridge 


Scour Data Management System (BSDMS) Coarse dataset. 


historical data from sites 60-62, which supports Mueller’s 
(1996) envelope curve for K, for lower velocity-intensity ratios 
(fig. 10). Nearly all data collected for this study plot below the 
idealized K,-envelope curve; one data point plots on the curve 
and six data points for five different sites (sites 20, 37, 51, 55, 
and 57) plot slightly above the curve. Idealized K,, values plot- 
ting above the envelope curve indicate that current methods 
(2011) underpredict pier scour and that factors not described 
by current (2011) equations led to scour underprediction. This 
underpredicted scour depth had a median residual, or differ- 
ence between measured and computed scour depth, of about 
0.2 ft. The most likely factor causing scour underprediction 
for the five sites where pier scour was measured for this study 
was the presence of woody debris on piers, noted as a moder- 
ate effect at four of the five sites (table 2). While this and other 
factors were not considered substantial enough to remove 
measurements from the dataset, the factors very likely con- 
tributed to the pier-scour depth measured in the field but not 
predicted by the HEC-18 equations (fig. 10). Additionally, five 
sites from the BSDMS Coarse dataset show underpredicted 
scour with a median residual of 0.41 ft (fig. 10). The five 
BSDMS sites were discussed by Mueller and Wagner (1995), 
where four of the five idealized K, values were found to plot 
below the envelope line when the D,, and D,. variables used in 
the prediction equations were derived from composite samples 
of bed-material obtained at each site. 

Mueller and Wagner (2005) proposed a different equation 
for computing K, (equation 20), which related K, and rela- 
tive bed-material size. When plotted, the equation defines an 


envelope curve that represents the maximum idealized K, for a 
given bed-material size. Idealized K,, and relative bed-material 
sizes for the pier-scour measurements collected for this study, 
measurements from the BSDMS Coarse dataset, and histori- 
cal pier-scour measurements from sites 60-62 were plotted 
along with the equation Mueller and Wagner (2005) proposed 
(fig. 11). All data plot below the K,-envelope curve proposed 
by Mueller and Wagner (2005), indicating that scour was not 
underpredicted. The lower end of the relative bed-material size 
scale is mostly populated with idealized K, values computed 
from data collected for this study because of the larger bed- 
material sizes involved (fig. 11). The relation between ideal- 
ized K, and relative bed-material size are similar for data col- 
lected for this study, the BSDMS Coarse dataset, and historical 
data for sites 60—62 and therefore, all data support the Mueller 
and Wagner (2005) envelope curve because pier scour was 
always overpredicted. 

The maximum value of relative pier scour (y’,/b) 
obtained in the study was 1.00 (table 5), and the maximum 
values observed in the two BSDMS datasets (All, Coarse) 
were 2.09 and 1.25, respectively. Research indicates that 
maximum scour generally does not exceed about 2.4 times 
the effective pier width (y’,= 2.45) for measurements with 
Froude numbers less than about 0.8 (Richardson and Davis, 
2001), which is supported by the data collected for this study 
and BSDMS datasets. The largest maximum of 2.09 (BSDMS 
All dataset) includes data for sites with sand bed material, 
whereas the lower maxima (1.00 and 1.25) include data col- 
lected for this study and from the BSDMS Coarse dataset for 
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5 T T TTT T TT TTTTTyT T TT TTI T T 1 TTTITy T T .TTTTTl EXPLANATION 


@ Data collected for this study (n=103) 


v Data from sites 60-62 for 1996-97 (n=9) associated with 
peak streamflows of greater recurrence interval than 
data collected for this study 


wn 
a Envelope curve for K, as proposed by X Data from BSDMS Coarse dataset (=90) 
> Mueller (1996) and applied in current 
=} te HEC-18 criteria (Richardson and L = ctl Variables 
g r Davis, 2001) a 4 
a L : _- % 4 n_ is the number of measurements 
= L es | 
a L x x aw 4 K, is the correction factor for the coarse bed-material 
= L @ x= 4 armoring effect on pier-scour depth, dimensionless 
SS x X@ ~ ox xX x 
a oe e Ox @ x 4 oo hae: . ; 
& ‘ X Qe® 3g x 7 x J, 1s pier-scour depth from field measurements, in feet 
os L Ge ®@ 9 $8. as x 56 | y, is pier-scour depth computed by an equation, in feet 
Ss L @ xe . eX pe x J V, is the approach velocity directly upstream from 
= oe. ori 
@ QR the pier, in feet per second 
a fo) e %O OPE @ s is 
% « e@ xx Ms9_is the critical velocity for incipient motion of bed material 
e ® i : aun 
e ‘ upstream of the bridge for the D., particle size, in feet 
01 - fo) xXx x ov 7 per second 

- 4 = 

F 4 Vso is the critical velocity to initiate scour of bed material at 

r AA 7 he pier for the D., particle size, in feet per second 

0.05 1 Popririiil i piiiriiit I pi iiiitt 1 pouririil 1 Pouriiit 
0.01 0.1 1 10 100 1,000 Vio is the critical velocity to initiate scour of bed material at 


he pier for the D,, particle size, in feet per second 
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| DIMENSIONLESS 
Figure 10. Comparison of idealized K, to Mueller velocity-intensity ratios for pier-scour measurements collected for this study, 


historical pier-scour measurements from Montana, and pier-scour measurements in the Bridge Scour Data Management System 
(BSDMS). 
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Data collected for this study (7=103) 


Data from sites 60-62 for 1996-97 (n=9) associated with peak 
4 streamflows of greater recurrence interval than data 


0.19 
Ky = 0.35 (b/Ds9) collected for this study 


Envelope curve proposed by 
Mueller and Wagner (2005) Xx Data from BSDMS Coarse dataset (n=90) 


Variables 


n is the number of measurements 


| Ky is the correction factor for the coarse bed-material armoring 
7 effect on pier-scour depth, dimensionless 


' 


IDEALIZED K, (’s /y, ), DIMENSIONLESS 


y' is pier-scour depth from field measurements, in feet 
¥, is the pier-scour depth computed by an equation, in feet 


J ob is effective width of the pier measured perpendicular to the 
approach flow, in feet 


Dsy _ is the bed-material particle size where 50 percent of 


0.1 particles are finer by weight, in millimeters or feet 


0.05 1 1 poi i ritt 1 1 1 
1 10 100 1,000 


RELATIVE BED-MATERIAL SIZE (b/Ds, ), DIMENSIONLESS 


Figure 11. Comparison of idealized K, to relative bed-material size for pier-scour measurements collected for this study, historical pier- 
scour measurements from Montana, and pier-scour measurements in the Bridge Scour Data Management System (BSDMS) Coarse 
dataset. 


sites characterized by much coarser bed material. These results 
reaffirm the general finding that pier scour is inversely related 
to bed-material size. 


Comparisons of Measured and Computed Scour 


The most reliable pier-scour equation minimizes under- 
predicted scour depth (positive residual) as the principal 
criterion but also does not unreasonably overpredict scour 
depth (negative residual). Scour depth was computed for the 
measurements collected in this study using the basic HEC-18 
equation and five pier-scour equations described previously 
that use the K, correction factor for the armoring effect of 
coarse-bed material: 


1. The basic HEC-18 equation with no K, factor (equa- 
tion 3 or 4). This is the original equation from which 
other equations have evolved and is the equation 
(equation 4) used in comparisons involving the ide- 
alized K, (equation 22). 


2. The HEC-18 equation (equation 5), with K, deter- 
mined by equation 11 as proposed by Mueller (1996, 
where 0.4 < K,< 1.0) and used by Richardson and 
Davis (2001). This equation is currently (2011) rec- 
ommended by the FHWA and is referred to hereinaf- 
ter as the HEC-18-K,Mu equation. 


3. The HEC-18-K,Mu equation (equation 5) with the 
current (2011) FHWA constraint (0.4 < K,< 1.0) 
removed and 0.0 < K,< 1.0 is applied instead. Values 
of K, determined by equation 11 to be less than 0.4 
are, therefore, used in equation 5. This version of 
equation 5 is referred to hereinafter as the HEC-18- 
K,Mu' equation. 


4. The equation developed using nonuniform sediment 
mixtures (equation 13) by Molinas (2003), referred 
to hereinafter as the Molinas (2003) equation. 


5. A variation of the HEC-18 equation (equation 19) 
proposed by Molinas (2003), with K, determined by 
equation 17 and referred to hereinafter as the HEC- 
18-K,Mo equation. 


6. The HEC-18 equation (equation 5), with K, deter- 
mined by equation 20 as proposed by Mueller and 
Wagner (2005) and referred to hereinafter as the 
HEC-18-K,MW equation. 


Computed scour depths were compared to measured 
scour depths using scatter plots (fig. 12) and boxplots of 
residuals (fig. 13) to evaluate the reliability of the equations to 
predict scour. Additionally, positive and negative residuals of 
scour were computed and reported here because of the impli- 
cation each type of residual separately has to bridge founda- 
tion design. The evaluation showed: 
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¢ The basic HEC-18 equation (fig. 124 and fig. 13) did 
not underpredict scour depth at any site, but scour 
depth was substantially overpredicted (median and 
maximum overpredicted residuals of about -4.0 ft and 
-12.4 ft, respectively). 


The HEC-18-K,Mu equation (fig. 12 and fig. 13) 
underpredicted scour depth six times at five different 
sites (median and maximum underpredicted residuals 
of about 0.2 ft and 0.5 ft, respectively) and moderately 
overpredicted scour depth (median and maximum 
overpredicted residuals of about -0.9 ft and -4.4 ft, 
respectively). 


The HEC-18-K,Mu’ equation (fig. 12C and fig. 13) 
underpredicted scour depth 12 times at 9 different sites 
(median and maximum underpredicted residuals of 
about 0.3 ft and 0.8 ft, respectively) and moderately 
overpredicted scour depth (median and maximum 
overpredicted residuals of about -0.7 ft and -4.4 ft, 
respectively) similar to the HEC-18-K,Mu equation. 


The Molinas (2003) equation (fig. 12D and fig. 13) 
underpredicted scour depth 35 times at 23 different 
sites (median and maximum underpredicted residuals 
of about 0.4 ft and 2.0 ft, respectively) and moderately 
overpredicted scour depth (median and maximum 
overpredicted residuals of about -0.6 ft and -2.8 ft, 
respectively) similar to the HEC-18-K,Mu and HEC- 
18-K,Mu’ equations. 


The HEC-18-K,Mo equation (fig. 12£ and fig. 13) 
slightly underpredicted scour depth at three sites 
(median and maximum underpredicted residuals of 
about 0.1 ft and 0.2 ft, respectively), but overpre- 
dicted scour depth more substantially (median and 
maximum overpredicted residuals of about -2.0 ft and 
-8.6 ft, respectively) compared to the HEC-18-K,Mu, 
HEC-18-K,Mu’, and Molinas (2003) equations. 


The HEC-18-K,MW equation (fig. 12/ and fig. 13) 
did not underpredict scour depth at any site, but scour 
depth was overpredicted (median and maximum over- 
predicted residuals of about -1.9 ft and -9.9 ft, respec- 
tively) similar to the HEC-18-K,Mo equation. 


Finally, boxplots summarizing residuals for the six equa- 
tions (fig. 13) are consistent with the conclusion of Chase and 
Holnbeck (2004), which indicated that the HEC-18-K,Mu 
equation generally predicted scour depth in closer agreement 
to measured scour than the other equations used for computing 
pier-scour depth in coarse-bed streams. Furthermore, the mag- 
nitude of the residuals of underpredicted and overpredicted 
scour depth was the lowest using the HEC-18-K,Mu equation. 
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were not plotted | 
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were not plotted 


© HEC-18-K,Mo 


© HEC-18-K,MW 


COMPUTED PIER-SCOUR DEPTH (y,), IN FEET 
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0 0 1 1 1 i 1 1 1 i 1 1 1 i 1 1 1 i 1 1 1 
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MEASURED PIER-SCOUR DEPTH (y’,), IN FEET MEASURED PIER-SCOUR DEPTH (y’,), IN FEET 
EXPLANATION 
<_ isless than y, is the pier-scour depth computed by an equation, in feet 
> is greater than y’, is pier-scour depth from field measurements, in feet 
< 


is less than or equal to K 


, isthe correction factor for the coarse bed-material armoring 


effect on pier-scour depth, dimensionless 


Figure 12. Comparison of measured with computed pier-scour depths for data collected for this study using selected pier-scour 
equations. A, Using the basic HEC-18 equation without K,. B, Using the HEC-18-K,Mu equation, where 0.4 < K, < 1.0. C, Using the 
HEC-18-&,Mu’ equation, where 0.0<K, < 1.0. D, Using the Molinas (2003) equation. £, Using the HEC-18-,Mo equation. 

F, Using the HEC-18-K, MW equation. 
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ing effect on pier-scour 
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Residual pier-scour depth is measured 
pier-scour depth minus computed 
pier-scour depth. 


Figure 13. Comparison of residual pier-scour depth for six pier-scour equations evaluated for this study. 


Effects of Bed-Material Characteristics on Pier 
Scour 


Bed-material data were collected for the surface or armor 
layer as part of each measurement. In addition, similar data 
also were collected for bed material in the shallow-subsurface 
bed material immediately below the surface layer. The 
shallow-subsurface layer included only that part of the stream- 
bed where scour and depositional processes are active (active 
layer). These detailed bed-material data for each layer (table 9) 
were used to assess the effects of bed-material characteristics 
on pier scour in coarse-bed streams in Montana. 

Assessment of the effects of bed-material characteristics 
on pier scour began with a detailed examination of the rela- 
tion of bed-material size and gradation of the surface layer to 
relative pier scour. This effort largely involved applying the 
work of previous researchers to the pier-scour measurements 
made on coarse-bed streams for this study. Next, the effects of 
shallow-subsurface bed-material characteristics on pier scour 
were assessed and compared to the effects of the surface bed- 
material characteristics. In particular, the differences in pier- 
scour prediction caused by using bed-material characteristics 
of one layer or the other were examined. Finally, the effects 
of differences in bed-material size and gradation between the 
layers on pier scour were each assessed separately and then 
the combined effects of differences in bed-material size and 
gradation between the two layers were assessed. 


Effects of Surface Bed-Material Characteristics 


Early bridge scour studies did not consider surface 
bed-material characteristics, but in the mid—1960s research- 
ers began to note that surface bed-material size had an effect 
on pier-scour depth (Shen and others, 1966). Later studies 
showed that the gradation of the surface bed material also 
affected pier-scour depth (Ettema 1976, 1980; Chiew, 1984; 
Baker, 1986; and Abdou, 1993). The effect of bed-material 
size and gradation on scour depth was investigated using 
data collected for this study, the BSDMS All dataset, and the 
BSDMS Coarse dataset to confirm and extend the findings of 
previous studies. Data from these three datasets were com- 
pared to previously published relations of relative pier scour 
(y’,/b) and relative velocity (V,/V_.,). Relative velocity was 
computed based on the critical velocity for incipient motion 
of bed material upstream from the bridge (V,.,). Comparisons 
also were made using relative velocity computed based on 
the critical velocity to initiate scour of bed material at the pier 
(V"..,). This was considered appropriate because most pier- 
scour data collected for this study reflect clear-water scour 
conditions, where incipient motion of the bed material and 
substantial sediment supply to the scour hole was not a factor. 

Recognizing the relation between bed-material charac- 
teristics and pier-scour depth, researchers have continued to 
improve scour-prediction equations by incorporating math- 
ematical expressions that attempt to account for bed-material 


40 Investigation of Pier Scour in Coarse-Bed Streams in Montana, 2001 through 2007 


Table 9. Particle-size distribution and gradation coefficient for surface and shallow-subsurface bed-material data collected for this 
study. 


[All pier-scour measurement data were obtained for bed material that was noncohesive. Symbols: D,, bed-material particle size where i percent of particles are 
finer by weight, in millimeters; 6.4 gradation coefficient of the surface bed material, dimensionless; Ons gradation coefficient of the shallow-subsurface bed 
material, dimensionless] 


Site Surface bed material 
number Site name 
(fig.1) D,, D,, D., D,, Ds Dy D,, D,, Gs 
Missouri River Basin 
1 Big Hole River east 40.5 73.3 102 176 179 209 263 309 2.09 
of Wise River 
2 Big Hole River west DBT 51.9 les ISS) 159 177 190 200 D5 
of Divide 
3 Beaverhead River 20.3 28.8 35.8 59.2 60.1 66.2 80.3 103 1.71 
north of Dillon 
4 Big Hole River near 28.2 42.7 SehS 104 107 123 157 Wo 1892; 
Melrose 
5 Boulder River near 29.1 44.4 58.0 111 114 127 163 284 1.96 
Basin 
6 Wisconsin Creek near DB) 36.0 47.8 78.8 80.2 SuEe 104 123 1.84 
Sheridan 
7 Gallatin River 32.7 51.2 64.7 117 118 126 167 272 1.89 
northwest of West 
Yellowstone 
8 Big Hole River 20.8 30.6 38.5 60.9 61.6 68.2 80.1 89.6 ail 
southwest of Twin 
Bridges 
9 Jefferson River north 28.8 44.9 60.7 114 116 125 148 174 1.99 
of Twin Bridges 
10 Madison River south 29.9 56.7 79.8 149 155 206 253 291 DIB 
of Cameron 
11 Boulder River at I-90, 2.58 10.3 17.1 44.1 45.8 61.8 82.9 188 4.14 
near Cardwell 
12 Boulder River near 2.58 10.3 17.1 44.1 45.8 61.8 82.9 188 4.14 
Cardwell 
13 South Boulder River 29.3 43.2 55.2 111 114 129 155 175 1.95 
near Cardwell 
14 South Willow Creek 8.76 29.1 38.3 77.3 79.0 87.0 109 151 2.97 
near Harrison 
15 Gallatin River west 45.5 58.5 70.4 109 112 127 143 179 1.55 
of Bozeman 
16 Jefferson River west 5.91 12.4 2.3} 57.0 59.3 73.6 89.6 126 3.11 
of Three Forks 
17 Gallatin River near 17.8 26.7 33.3 57.4 58.3 63.9 76.0 102 1.80 
Logan 
18 Jefferson River north 8.85 16.1 19.1 31.1 31.5 38.2 47.5 63.7 1.87 
of Three Forks 
19 Madison River near 19.7 25.3 29.8 43.2 43.7 48.2 56.8 63.7 1.48 
Three Forks 
20 Missouri River near 11.9 19.9 29.7 64.1 66.2 76.8 87.4 118 D3) 
Townsend 
21 Sun River north of 26.8 45.7 57.4 111 114 128 156 179 2.03 


Augusta 
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Shallow-subsurface bed material 


D D. D. D. D D. D. D. 6 


16 35 50 84 85 90 95 99 GSS 


Missouri River Basin 


11.8 39.1 66.4 139 142 154 167 177 3.44 
V9) Sil Sts} 109 110 116 ee 127 3.86 
3.26 12.9 21.8 51.8 53.2 60.4 72.6 86.5 3.99 
4.95 18.7 47.5 104 106 113 121 127 4.59 
4.35 22.4 40.4 74.4 75.4 80.3 85.1 89.0 4.14 
1.97 17.3 31.1 71.8 73.0 78.6 84.3 88.9 6.03 
6.06 24.5 41.3 73.8 74.9 79.9 85.0 89.0 3.49 
6.43 19.6 DOs) 48.5 49.4 54.3 Se). 63.0 Dele) 

14.3 38.9 yet) 86.6 87.4 95.7 112 125 2.46 

Bei 93.8 119 161 163 168 174 179 2.66 
1.03 1.66 3.64 27.6 28.5 33.5 41.4 57.0 5.18 
1.03 1.66 3.64 DKS BSS) 3)3).5) 41.4 57.0 Sats 
6.89 29.5 53.8 108 109 116 122 127 3.96 
1.87 2,1 32.8 105 106 114 121 127 7.49 
6.05 28.1 58.4 147 149 159 170 178 4.93 
Bell 10.4 16.0 36.6 38.0 45.0 54.5 62.1 3.64 
9.24 25.2 35.4 67.4 68.9 73:9 83.0 88.6 2.70 
1.86 5.44 O19? 24.9 BSS) 30.9 40.1 56.9 3.66 
1.91 14.5 19.0 35.0 36.3 42.8 52.6 61.7 4.28 
1.85 0.38) 18.2 40.9 41.7 45.6 54.8 62.2 4.71 


2.62 12.4 32.4 153 154 163 171 178 7.63 
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Table 9. Particle-size distribution and gradation coefficient for surface and shallow-subsurface bed-material data collected for this 
study.—Continued 


[All pier-scour measurement data were obtained for bed material that was noncohesive. Symbols: D,, bed-material particle size where i percent of particles are 
finer by weight, in millimeters; 6, , gradation coefficient of the surface bed material, dimensionless; Ons gradation coefficient of the shallow-subsurface bed 
material, dimensionless] 


Site Surface bed material 
number Site name 
(fig.1) D,, D,, D., D,, Dy Dy D,, D,, Gis 
Missouri River Basin—Continued 

22 Marias River near 7.36 17.2 22.7 44.0 45.0 52.1 59.3 T4.4 2.44 
Shelby 

2B Smith River southeast 6.67 10.4 13.6 26.7 DUS 31.2 41.5 58.7 2.00 
of Ulm 

24 Smith River at old 19.8 30.9 39.2 86.0 89.7 106 123 165 2.08 


truss bridge, south- 
east of Ulm 


De) Belt Creek 7 miles Sel 29.6 44.6 88.8 90.2 108 125 168 2.42 
south of Belt 

26 Belt Creek 6 miles 14.5 28.1 38.0 75.1 77.0 86.5 105 123 2.28 
south of Belt 

27 Belt Creek 2 miles MBM 3133 39.5 76.3 oe, 87.9 109 OF 1.82 


southeast of Belt 


Yellowstone River Basin 


28 Yellowstone River at 14.6 25.0 3352 67.6 69.3 78.9 94.4 110 DMS 
Emigrant 

29 Shields River near 24.8 40.2 50.6 88.0 89.6 102 115 125 1.88 
Wilsall 

30 Mill Creek near Pray! 32.2 Sil 67.4 130 133 149 164 177 2.01 

31 Yellowstone River 41.5 71.5 92.1 133 136 150 165 177 1.79 
near Pray! 

32 Yellowstone River 10.9 330) 2) 54.5 114 116 126 174 244 3328} 
near Pine Creek 

33 Yellowstone River 24.1 44.5 55.5 88.8 89.9 103 117 128 1.92 
south of Livingston 

34 Boulder River near 31.9 42.9 55.6 108 110 119 128 178 1.84 
McLeod 

35 Boulder River south 43.2 61.3 76.6 125 127 148 171 269 1.71 
of Big Timber 

36 Boulder River at I-90 46.9 67.7 84.6 136 140 164 196 266 1.70 


near Big Timber 


Columbia River Basin 


37 Little Blackfoot River ie) 23 29.2 52} 56.4 62.6 BS) 119 2.06 
near Avon 


38 Blackfoot River (old 1.68 4.31 9.63 26.8 27.6 31.6 38.2 43.6 4.00 
bridge) west of 
Lincoln! 


39 Blackfoot River (new 1.68 4.31 9.63 26.8 27.6 31.6 38.2 43.6 4.00 
bridge) west of 
Lincoln! 


40 Little Blackfoot River 14.6 42.6 57.4 98.1 101 114 127 169 2.59 
near Garrison 
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Shallow-subsurface bed material 


D, D,, D., D,, D,, D,, D,, Dg o. 
Missouri River Basin—Continued 
4.54 14.9 25.3 65.2 66.8 74.5 82.3 88.5 3.19 
1.59 3.0) 9.10 30.6 31.6 378) 44.3 59.8 4.39 
10.4 20.1 29.0 55.0 56.0 60.8 70.4 86.1 2.30 
1.98 Weil 18.0 63.8 65.3 73.6 81.8 88.4 5.68 
1.70 14.6 27.2 78.5 79.2 82.8 86.4 89.3 6.79 
2.00 9.41 20.0 67.3 69.0 75.8 82.9 88.6 5.80 
Yellowstone River Basin 
1.94 15.3 ATS) 60.4 61.2 67.7 78.9 87.8 S57) 
2.30 10.3 19.3 68.8 70.3 78.7 87.3 94.3 5.47 
3) 51.0 67.4 130 133 149 164 WaT 2.01 
41.5 72.0 92.1 133 136 150 165 177 1.79 
5.74 2S S91 151 152 162 171 178 5.13 
7.06 22.5 42.2 147 149 159 170 178 4.56 
8.86 48.5 67.8 110 111 117 122 127 Brod) 
7A1 38.0 66.4 104 106 113 121 127 3.83 
13.8 56.8 101 158 159 166 173 179 3.38 
Columbia River Basin 
1.99 6.80 14.8 46.0 47.5 54.9 62.3 83.3 4.81 
1.68 4.31 9.63 26.8 27.6 31.6 38.2 43.6 4.00 
1.68 431 9.63 26.8 27.6 31.6 38.2 43.6 4.00 
9.39 20.1 31.3 58.1 58.8 62.4 73.3 86.7 2.49 
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Table 9. Particle-size distribution and gradation coefficient for surface and shallow-subsurface bed-material data collected for this 
study.—Continued 


[All pier-scour measurement data were obtained for bed material that was noncohesive. Symbols: D,, bed-material particle size where i percent of particles are 
finer by weight, in millimeters; C4 gradation coefficient of the surface bed material, dimensionless; Ons gradation coefficient of the shallow-subsurface bed 
material, dimensionless] 


Site Surface bed material 
number Site name 
(fig.1) D,, D,, D,, D,, De Dy D,, D,, Gis 
Columbia River Basin—Continued 

41 Clark Fork near Gold 8.62 19.8 30.9 79.8 82.1 97.2 159 165 3.04 
Creek 

42 North Fork Black- 22.4 29.8 36.6 60.6 61.5 69.7 82.2 108 1.65 
foot River west of 
Lincoln 

43 North Fork Blackfoot 25.9 39.6 50.1 93.7 96.8 112 128 169 1.90 
River near Ovando 

44 Clark Fork south of 5.46 14.4 34.4 86.4 87.9 102 118 153 3.98 
Drummond 

45 Lost Horse Creek 24.1 52.0 73.0 122 123 139 168 264 2.26 
north of Darby 

46 Clark Fork near 15.8 21.8 27.8 43.6 44.0 51.4 62.0 83.6 1.66 
Drummond 

47 Clark Fork at Turah 20.9 37.1 49.3 87.1 88.3 103 122 165 2.04 
Bridge 

48 Swan River east of 34.8 RYO) 72.9 LS) 117 124 146 733 1.82 
Ferndale 

49 Clark Fork at Mis- 28.5 47.2 67.3 123 125 144 168 197 2.08 
soula 

50 Clark Fork near 18.1 29.7 46.1 95.4 98.9 116 140 ee 2.30 
Alberton 

51 Clark Fork near 20.0 36.9 53.4 105 108 121 144 174 2.29 
Superior 

Sy) St. Regis River at Dor 35.5 50.7 88.6 89.8 105 121 161 1.97 
Ward Creek 

53 St. Regis River west 27.2 35.8 41.0 69.6 71.1 78.8 86.4 114 1.60 
of St. Regis 

54 St. Regis River at St. 17 19.5 24.5 50.0 SP) 63.2 84.8 OF 1.99 
Regis 

55 Flathead River near 16.4 21.4 27.6 54.7 56.3 64.3 81.6 125 1.82 
Perma 

56 Clark Fork above 16.2 38.4 54.9 115 117 128 159 240 2.67 
Flathead River, 
near Paradise 

57 Clark Fork below 18.7 255 31.3 51.1 52.4 58.8 69.4 102 1.65 
Flathead River, 
near Paradise 

58 Kootenai River near 27.4 45.8 58.0 93.0 96.0 111 126 168 1.84 
Libby 

59 Yaak River near Troy 46.1 80.6 108 189 196 231 265 293 2.02 


'The particle size of the shallow-subsurface bed material was estimated in the field to be essentially the same as the surface bed material. 
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Shallow-subsurface bed material 


D,, D,, D., D,, D,, D,, D,, Dy, o. 
Columbia River Basin—Continued 
1:92 17.5 31:5 72.2 73.3 78.8 84.4 88.9 6.13 
1.49 5.08 13.2 44.0 45.1 51.4 Sod 62.7 5.44 
3.30 10.7 18.5 61.9 63.9 72.6 81.3 88.3 4.33 
4.25 23.4 40.5 76.4 WIS 81.5 85.8 89.2 4.24 
11.4 46.4 69.1 100 102 111 119 126 2.96 
4.69 13.2 19.3 45.9 47.0 52,7) 58.3 62.9 33 
3.20 20.8 33.7 62.2 62.9 70.8 80.4 88.1 4.41 
8.60 30.6 S27) 107 108 115 121 127 3),.5) 
4.09 18.9 36.0 136 138 152 166 177 5.76 
2.40 10.7 21.8 76.8 78.8 88.8 108 124 5.66 
5.31 21.0 52.3 97.1 99.1 109 118 126 4.27 
25.0 62.4 Qi Il 116 aw 121 124 27 2.16 
1.90 7.20 15.2 36.9 37.4 40.0 42.5 44.5 4.41 
2.49 16.3 320) 88.0 89.7 102 115 125 5.93} 
17 16.4 26.9 72.1 73.2 78.8 84.4 88.9 3.73 
5.95 20.8 39.1 77.6 78.3 82.2 86.1 89.2 3.61 
1.78 13.0 17.4 38.5 39.5 44.3 76.0 87.2 4.65 
12.1 36.0 66.4 86.8 87.4 92.4 110 124 2.68 
17.9 47.9 73.1 111 112 117 123 127 2.49 
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effects, as cited in the section “Development of Pier-Scour 
Equations for Coarse-Bed Streams.” Some of these improve- 
ments were based on laboratory studies and some were based 
on field studies. Many of the equations developed or improved 
based on laboratory studies are not able to successfully predict 
pier-scour depth in natural streams with highly variable bed- 
material characteristics. Even equations developed based on 
field data are not able to successfully predict pier scour under 
all field conditions. For example, the most recent (2011) 
HEC-18 equation (equation 5) incorporates a correction factor 
for coarse bed streams (K,). K, is calculated using a velocity- 
intensity ratio (equation 10) that is dependent on the D,, and 
D,, particle sizes. Under most natural conditions, the D,. 
moderately exceeds the D.,, the velocity-intensity ratio and 

K, can be calculated, and the resulting scour predictions are 
acceptable. Mueller and Jones (1999, p. 306) and Richardson 
and Davis (2001, p. 6.6) reported that the configuration of the 
denominator in equation 10 can produce anomalous results 
(negative velocity ratio) for rare combinations of D,, and D,.. 
None of the surface bed-material particle-size data associated 
with pier-scour measurements collected for this study pro- 
duced anomalous results when used in equation 10. To better 
understand conditions where equation 10 produced anomalous 
results, selected pier-scour data collected for this study were 
used in equation 10 and manipulated (increased or decreased) 
until anomalous results were achieved. Anomalous results 
were produced when D,. substantially exceeded D,, and when 
D,, was similar in magnitude to D,, for certain D,, and D,. size 
combinations. Because the solution of equation 10 depends on 
D,, being larger in some modest proportion to D,,, the surface 
bed-material data from this study were analyzed to statistically 
determine the relation between the D,, and D,. particle sizes 
based on regression analysis (Helsel and Hirsch, 2002) and 
evaluation of r? and p. The coefficient of determination (r’) 
has a value that ranges from 0 to | and the strength of the rela- 
tion increases as the value of r* approaches 1. The p-value is 
the probability of obtaining a result that is more extreme than 
a given value, where a p-value less than 0.05 indicates statisti- 
cal significance (significance level) of the regression relation. 
Based on the regression analysis, the D,, and D,. particle sizes 
have a positive relation with each other that is statistically sig- 
nificant (r* = 0.73 and p < 0.001). In other words, for the data 
collected for this study, sites that have a larger D,, particle size 
also tend to have a larger D,. particle size. This observation is 
similar to the observation of Mueller and Jones (1999) about 
the bed-material data originally used to develop the equation 
for the K, correction factor. 

Bed-material size is an important characteristic that deter- 
mines the magnitude of the velocity required to cause pier 
scour. The effect of bed-material size on pier scour typically 
is investigated in conjunction with hydraulic characteristics 
like velocity and flow depth because these characteristics 
also influence the mobility of bed material. Relative velocity 
(V/V ,,) incorporates all three of these characteristics. Previ- 
ous researchers have examined the relation between relative 
velocity and relative pier scour and formulated envelope 


curves. Melville (1984) and Chiew (1984) found the first peak 
in the envelope curve to be on the threshold between clear- 
water and live-bed scour (V,/V_,, = 1.00) and a second live-bed 
scour peak became asymptotic (approaching a constant value) 
with increasing relative velocity and live-bed scour condi- 
tions. They also showed that the magnitude of the second peak 
was less than the first peak for nonripple-forming sediments 
(D,,2 0.6 mm), and the second peak was greater than the first 
peak for ripple-forming sediments (D,,< 0.6 mm). 

The relation between relative pier scour and relative 
velocity defined by an envelope curve for the data collected 
for this study and the BSDMS Coarse dataset (fig. 144) were 
compared to relations observed in laboratory data by Mel- 
ville (1984) and Chiew (1984). The envelope curve for data 
used for this study displays an earlier peak with a decrease 
in relative pier scour followed by an asymptotic rise. Inter- 
estingly, the first peak (fig. 144) in the curve occurs at a 
lower relative velocity (V./V_,, ~ 0.75) than the clear-water 
and live-bed threshold (V,/V_., = 1.00). Also, the magnitude 
of relative pier scour (maximum of 1.25) is much less than 
the published value of 2.4 obtained from laboratory data for 
cylindrical piers (Melville and Sutherland, 1988). Further- 
more, equilibrium live-bed scour is about 18 percent shallower 
than the maximum pier scour under clear-water conditions, 

a larger reduction than the 10 percent indicated by Richard- 
son and Davis (1995, p. 16). Montana historical pier-scour 
data from 1996-97 (fig. 144) plot below or on the envelope 
curve, indicating that the curve developed from scour data for 
flows approximating bankfull discharge (< 2-year recurrence 
interval) also envelopes data for more severe hydrologic and 
hydraulic conditions with recurrence intervals up to about 
100 years (table 8). 

As noted, the envelope curves that bound the pier-scour 
data collected for this study (fig. 144), and in particular the 
magnitude of the maximum relative pier scour and the relative 
velocity associated with that maximum, were different than 
expected based on envelope-curve relations from previous 
studies. Some of these differences between data and relations 
from past studies and data and relations from this study could 
be mostly reflections of the paucity of data defining the y’ /b 
maxima (fig. 144), with the true maxima being slightly greater 
and corresponding relative velocity (V,/V_.,) closer to the 
theoretical clear-water and live-bed scour threshold. However, 
differences in the bed-material characteristics and flow condi- 
tions between data collected for this study and earlier studies 
are a more likely explanation. First, the sizes of bed-material 
associated with pier-scour data collected for this study and 
included in the BSDMS Coarse dataset are much larger 
than bed-material sizes associated with pier-scour data from 
research (Melville and Sutherland, 1988) that documented 
higher maximum relative pier scour values approaching 2.4. 
Less pier scour would be expected at sites with coarser bed 
material and associated armoring. In addition, hydraulic 
conditions encountered in the field that cause scour of bed 
material at the pier typically involve variable streamflows over 
some duration (unsteady flow), compared to constant flow 
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EXPLANATION 
@ Data collected for this study (m=103) 


A Data from the Bridge Scour Data Management 
System (BSDMS) Coarse dataset (n=90) 


BH Data for sites 60-62 for 1996-97 (n=9) associated 
with peak streamflows of greater recurrence 
interval than data collected for this study 


Variables 
n__ is the number of measurements 


y's is pier-scour depth from field measurements, 
in feet 


b isthe effective width of the pier measured 
perpendicular to the approach flow, in feet 


Vo_ isthe approach velocity directly upstream from 
the pier, in feet per second 


is the critical velocity for incipient motion of bed 
material upstream of the bridge for the 
Ds0¢ particle size, in feet per second 


V'e50 is the critical velocity to initiate scour of bed 
material at the pier for the Dsq particle size, 
in feet per second 
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Figure 14. Relation between relative pier scour and: A, relative velocity based on critical velocity for incipient motion of bed 
material upstream from the bridge; and B, relative velocity based on critical velocity to initiate scour of bed material at the pier. 


conditions over a duration (steady flow) typically used in the 
laboratory. Raudkivi and Sutherland (1981) investigated the 
effects that steady and unsteady flow conditions and flow dura- 
tion had on the relation between pier scour and relative veloc- 
ity in the laboratory. They found that pier scour attained under 
unsteady-flow conditions was only about 50 to 70 percent 

of pier scour attained under steady-flow conditions. Mueller 
(1996) made the observation based on the work of Raudkivi 
and Sutherland (1981) that the nearer the relative velocity of 
the flood peak is to the threshold for live-bed conditions and 
the longer the duration of the flood wave, the nearer the scour 
depth approaches the steady-state scour depth corresponding 
to the peak discharge. The magnitude of relative pier scour 


values for data collected for this study and from the BSDMS 
Coarse dataset (fig. 144) reflects unsteady-flow conditions 
with peak flows of limited duration and is, therefore, consis- 
tent with findings of past research (Raudkivi and Sutherland, 
1981; Mueller, 1996, p. 70-72). Armoring associated with 
nonuniform coarse-bed material and unsteady-flow conditions 
in the field probably account for differences between envelope 
curves developed for this study and those from previous stud- 
ies conducted in the laboratory using fine-grained material and 
steady-flow conditions. 

The very coarse nature of the bed material associated 
with pier-scour data collected for this study and in the BSDMS 
Coarse dataset most likely also contributed to the narrow 
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range of relative velocity (V/V,., ) ratios (about 0.55 to 1.55, 
with one value near 2.4, fig. 144). Finer sediment has been 
shown (Chiew, 1984) to produce a wider range of V/V. 
ratios (about 0.5 to 5.0). The prevalence of clear-water condi- 
tions associated with pier-scour measurements collected for 
this study indicates that the critical velocity to initiate scour 
of bed material at the pier (V’_.,, equation 2) might be used in 
place of the critical velocity for incipient motion of bed mate- 
rial upstream from the bridge (V.,,) for demonstrating rela- 
tions between relative pier scour and relative velocity. Replac- 
ing V/V. with V/V" .,as a measure of velocity intensity for 
pier scour (fig. 148) produces similar results (fig. 144), but in 
addition, shows that maximum relative pier scour (as mea- 
sured by y’,/b) occurs when the approach velocity (V,) is about 
1.3 to 1.4 times the critical velocity to initiate scour of bed 
material at the pier (V’,.,). If the approach velocity increases 
beyond that point, relative pier scour decreases until approach 
velocity reaches a minimum of about 1.8 to 1.9 times V"..., 
then increases asymptotically to reach equilibrium scour under 
live-bed conditions. 

The threshold for incipient motion of the streambed that 
defines the boundary between clear-water and live-bed scour 
conditions (V,/V_.,= 1.00) in figure 148 was determined by 
regression analysis of corresponding pairs of V/V... and 
V/V" ., ratios computed from the pier-scour data collected 
for this study and the BSDMS Coarse dataset. The resulting 
regression equation was used to compute the ratio of V/V". 
that statistically corresponded to V/V. = 1.00, the boundary 
between clear-water and live-bed scour conditions. Based on 


data collected for this study and the BSDMS Coarse dataset, 
the relative velocity associated with the threshold between 
clear-water and live-bed scour conditions is 1.94 times the 
relative velocity associated with initiation of scour at the pier 
(V/V). 

In addition to surface bed-material size, previous stud- 
ies showed that surface bed-material gradation also affected 
pier-scour depth. According to Ettema (1976), when gradation 
coefficients (o, ) are more than about 1.4 to 1.5, scour holes 
can become armored limiting further pier scour. Gradation 
coefficients associated with surface bed-material data collected 
for this study ranged from 1.48 to 4.14, with a median of 2.01 
(table 5). Bed-material gradation associated with pier-scour 
measurements collected for this study was less variable than 
bed-material gradation associated with pier-scour data in the 
BSDMS Coarse dataset (fig. 15), and bed-material gradation 
for all the pier-scour measurements collected for this study 
was sufficient for streambed armoring. 

The effect of surface bed-material gradation on pier 
scour was assessed as part of this study by relating the grada- 
tion coefficient (o_ ) to relative pier scour (fig. 16), similar to 
Mueller and Wagner (2005). Ripple-forming (fine) sediments 
(D,, < 0.6 mm) and nonripple-forming (coarser) sediments 
(D,, 2 0.6 mm) each have different relations, but the maximum 
relative pier scour for both kinds of sediments is concentrated 
near a value of 6, between about 2.0 and 2.5. After reaching a 
maximum value, relative pier scour decreases as oO, increases. 
Relative pier scour is greater for ripple-forming sediments 
when the gradation coefficient is less than about 2.5, except 
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Figure 15. Comparison of gradation coefficient for surface bed material for pier-scour measurements collected for this 


study with the Bridge Scour Data Management System (BSDMS) Coarse dataset. 


perhaps for small gradation coefficients (less than about 

1.3) where the envelope curves are poorly defined. Relative 
pier-scour for ripple-forming sediments reaches a maximum 
when the gradation coefficient is slightly less than or about the 
same magnitude as for nonripple-forming sediments, but the 
envelope curve for ripple-forming sediments (n=88) declines 
more sharply with increasing gradation coefficient compared 
to the envelope curve for the nonripple-forming sediments 
(n=184). Maximum relative pier-scour is substantially less 

for nonripple-forming sediments from this study because 

the bed material was much coarser compared to the BSDMS 
nonripple-forming sediment data. The envelope-curve decline 
rate associated with data from this study is not compared to 
the BSDMS envelope-curve decline rates because data from 
this study having gradation coefficients greater than about 4.0 
are lacking. Mueller and Wagner (2005) noted that the larger 
reduction in relative pier scour for ripple-forming sediments 
having more gradation may be due to armoring caused by the 
coarser sediment fractions. The average median particle size 
(D,,) for data points defining each of the three curves (fig. 16) 
also indicates that maximum relative pier scour decreases with 
an increase in bed-material size. 

Gradation effects also were investigated by using o, as a 
grouping variable when relating relative pier scour to relative 
velocity. Scour data were grouped into selected o, intervals, 
and the average o, value within an interval was used to define 
each group. No obvious relation was observed between rela- 
tive pier scour and relative velocity for any of the oO, intervals 
when the data were plotted, which is consistent with the 
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conclusion of Mueller and Wagner (2005). However, Ettema 
(1980) showed that relative pier scour for 0, > 2.5 is reduced 
for both ripple- and nonripple-forming sediments. Based on 
that observation, data collected for this study and the BSDMS 
Coarse dataset were grouped by 0., and two envelope curves 
were hand drawn (fig. 17A) that distinguish between Nes 2.5 
and cy 2.5. As shown by Ettema (1980) for data collected for 
this study and the BSDMS Coarse dataset, maximum relative 
pier scour is greater when o, < 2.5. 

Relations between relative pier scour, relative velocity, 
and gradation coefficient also were investigated using the 
critical velocity to initiate scour of bed material at the pier 
(VY...) in place of the critical velocity for incipient motion of 
bed material upstream from the bridge (V_.,) to define rela- 
tive velocity (fig. 17B). With relative velocity expressed as 
V/V" ..,, differences in maximum relative pier scour between 
the two envelope curves (fig. 178) that distinguish between 
oe 2.5 and a 2.5 are the same as before (fig. 17A), but the 
maxima demonstrate better symmetry with each other and are 
located at about 1.2 to 1.4 times the relative velocity defined 
by V/V’... The threshold for incipient motion of bed material 
upstream from the bridge that defines the boundary between 
clear-water and live-bed scour conditions for each envelope 
curve in figure 17B was determined by regression using cor- 
responding ratios of V./V,., and V/V’... computed from pier- 
scour measurements collected for this study and the BSDMS 
Coarse dataset. The resulting regression equation was used to 
compute the ratio of V/V’. that statistically corresponded to 
V/V, = 1.00 for each envelope curve. Based on data from 
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Figure 16. Relation between relative pier scour and gradation coefficient for pier-scour measurements collected for this study, 
historical pier-scour data from Montana, and pier-scour measurements in the Bridge Scour Data Management System (BSDMS) All 


dataset in ripple-forming (D,, < 0.6 millimeter) and nonripple-forming (D, 


«9 = 0.6 millimeter) sediments. 
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Figure 17. Relation between relative pier scour and: A, relative velocity based on critical velocity for incipient motion of bed material 
upstream from the bridge; and B, relative velocity based on critical velocity to initiate scour of bed material at the pier. 


this study and the BSDMS Coarse dataset, the relative velocity 
associated with the threshold between clear-water and live-bed 
scour conditions is 1.88 times the relative velocity associated 
with initiation of scour at the pier when o, < 2.5 (r? = 0.98 and 
p< 0.001), and is 2.00 times the relative velocity associated 
with initiation of scour at the pier when o, > 2.5 (r? =0.99 

and p < 0.001). Coarse-bed material with a larger gradation 
coefficient (0, >2.5) attains live-bed conditions at a higher 
relative velocity, compared to similar coarse-bed material with 
a smaller gradation coefficient (o, < 2.5). This observation 

is consistent with previous studies (Chiew, 1984 and Baker, 
1986) that noted the tendency for well-graded bed material 

to armor at lower relative velocities and resist further scour. 
Similarities in the magnitude of maximum relative scour for 
the envelope curve (fig. 178) associated with larger grada- 
tion coefficients (o, > 2.5) also may indicate that distinctions 
between clear-water and live-bed scour may not be important 
when surface bed material is coarse and well graded (9, > 2.5). 

Comparison of the two envelope curves indicates that a 
relation between relative pier scour and relative velocity that 
depends on o, may exist under field conditions but less conclu- 
sively than was indicated in laboratory research (Baker, 1986), 
where D,, was held constant while the size of the coarser 
fraction of the sediments was increased. A lack of field data 
for ripple-forming sediments with large gradation coefficients 
indicates that surface bed material having this characteristic 
may not be commonly found in natural streams. 

Pier-scour measurements collected for this study and the 
BSDMS AIl dataset confirm a general lack of data for coarse- 
bed sites with O, greater than about 2.5 (figs. 15 and 16). 
Grouping these data by size classification based on the D,, of 
the bed material also illustrates the lack of data for coarse-bed 
sites with 0, greater than 2.5 (table 10). The percentage of 
scour data from sites (n=375) with coarse-bed material (larger 
than sand, table 10) having Oo, greater than 2.5 is only 8 per- 
cent of the data. Furthermore, only 15 percent of scour data 
from sites with bed material classified as sand (table 10) have 
gradation coefficients exceeding 2.5. Mueller and Wagner 
(2005) found that the small number of sites with ripple-form- 
ing sediments with larger gradation coefficients made it dif- 
ficult to draw conclusions on the effect of gradation on scour. 
Similarly, the relatively small amount of pier-scour data from 
sites with coarse bed material and large gradation coefficients 
makes it difficult to understand the effects of gradation under 
these conditions. The narrow range in gradation coefficients 
displayed by all available data and the nonuniform conditions 
present in the field continue to complicate efforts to isolate the 
effects of gradation on pier scour. 
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Table 10. Percentage of pier-scour 
measurements' with a gradation coefficient 
greater than 2.5. 


[Abbreviations: BSDMS, Bridge Scour Data Man- 
agement System. Symbols: O,5 gradation coefficient, 
dimensionless; >, greater than; K,, correction factor 
for bed-material armoring effect; n, number of data 
values in dataset] 


Description class for Percent of data? 


surface bed material? with 6, >2.5 
Small cobble 1 
Very coarse gravel 2 
Coarse gravel 3 
Medium gravel 2 
Fine to very fine gravel 0 
Sand 15 


‘Includes all pier-scour measurements collected 
for this study (n=103) and BSDMS AI dataset 
(n=272) used in comparisons that meet HEC-18 cri- 
teria (Richardson and Davis, 2001) for applying K,. 


*Surface bed material classified based on the D,, 
particle size. 


3Number of observations within a description class 
(with o,> 2.5) divided by the total number of data 
(n=375), expressed as a percent. 


Effects of Shallow-Subsurface Bed-Material 
Characteristics 


Since the mid—1960s, pier-scour researchers have consid- 
ered the effects of bed-material characteristics on pier scour 
as cited in the section “Development of Pier-Scour Equations 
for Coarse-Bed Streams.” That research focused primarily 
on surface bed material because vertical variations in bed- 
material characteristics were not a concern in flume studies 
where the experimental bed material was composed of either 
uniform particles or a nonuniform mixture of sediment with 
the same size and gradation characteristics throughout the 
thickness of the bed. While collecting coarse-bed pier-scour 
data in the field for this study, differences were observed in 
surface and shallow-subsurface bed-material size and grada- 
tion (fig. 4). As a result, bed-material samples for both the 
surface and shallow-subsurface layers were collected and 
analyzed (table 9). The particle-size distribution and gradation 
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for the bed-material samples for each layer were compared 
and the differences in particle size that might relate to armor- 
ing were described by using ratios of the particle sizes for 
each layer. Differences in particle size and gradation between 
the two layers were apparent (table 9). Therefore, the effect 
these differences might have on pier-scour computations was 
investigated. 

Particle-size distribution curves for paired samples of 
surface and shallow-subsurface bed material collected for this 
study were analyzed and the D,, and D,. particle sizes were 
compared. The difference between the surface and shallow- 
subsurface bed-material D,, and D,. particle sizes (fig. 18) con- 
firmed what was observed in the field—due to armoring, larger 
particles were more commonly found on the streambed surface 
than in the shallow subsurface (fig. 4). However, coarse 
particles still made up a substantial portion of the shallow- 
subsurface bed material at most sites. Median D,, particle sizes 
for 103 surface and shallow-subsurface bed-material samples 
were about 49 mm and 32 mm, respectively, and median D,, 
particle sizes for 103 surface and shallow-subsurface bed 
material samples were 122 mm and 86 mm, respectively 
(fig. 18, table 5). These differences may result in reduced pier 
scour because bed-material size has been shown in this (figs. 6 
and 11) and other studies cited in the section “Development of 
Pier-Scour Equations for Coarse-Bed Streams” to be inversely 
related to pier-scour depth. 


Differences in gradation between surface and shallow- 
subsurface bed-material samples also were examined for this 
study. The median gradation coefficients for surface and shal- 
low-subsurface bed material associated with measurements 
made for this study were 2.01 and 4.14, respectively (fig. 19 
and table 5). The median gradation coefficient of the shallow- 
subsurface bed material was about twice that of the surface 
bed material (fig. 19) owing to the presence of finer-grained 
sediment in the shallow-subsurface bed material. These dif- 
ferences in gradation coefficient may result in reduced pier 
scour because gradation coefficients have been shown in this 
(figs. 16 and 17) and other studies cited in the section “Devel- 
opment of Pier-Scour Equations for Coarse-Bed Streams” to 
be inversely related to pier-scour depth in coarse-bed streams 
like those where pier-scour measurements were collected for 
this study. The combined effect of the larger surface bed- 
material size, relative to shallow-subsurface bed-material 
size, and the wider range of gradation coefficients in shallow- 
subsurface bed material is typical for the pier-scour measure- 
ments made for this study and may contribute to reduced pier 
scour in these coarse-bed streams. Differences in particle size 
and gradation coefficient between the two bed-material layers 
are likely similar to the variability in particle size and grada- 
tion coefficient that might be observed between multiple sets 
of pebble counts conducted to quantify average surface bed- 
material characteristics. This variability has been recognized 
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Figure 18. Comparison of surface and shallow-subsurface bed-material particle size for pier-scour measurements collected 


for this study. 
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SURFACE BED-MATERIAL SIZE (D. ,) IN MILLIMETERS 


EXPLANATION 
O  D,, particle size 
A Dys particle size 
r is the coefficient of determination 


p-value is the significance level 


Variable 
n is the number of measurements 
Dj. is the surface bed-material particle size, D, where 
i percent of particles are finer by weight, in millimeters 
or feet 
Dg, ss) 18 the shallow-subsurface bed-material particle size, D, 


where i percent of particles are finer by weight, in 
millimeters or feet 


Figure 20. Relations between corresponding surface and shallow-subsurface bed-material particle size (D,, and D,,) for pier-scour 
measurements collected for this study. 
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as a potential source of uncertainty in accurate prediction of 
pier-scour depth (Mueller and Wagner, 2005, p. 61). 

The relation between selected particle sizes of the two 
layers, independent of site location or specific particle size, 
was investigated using regression analysis of paired values 
of particle sizes for the surface and shallow-subsurface bed- 
material data collected for this study (fig. 20). A significant 
relation was indicated between surface and shallow-subsurface 
bed-material size for the D,, and D,. particle sizes (for D,,, 
r= 0.690 and p < 0.001; for D,., r? = 0.647 and p < 0.001). 
In other words, as the particles in the surface bed material 
increase in size, the particles in the shallow-subsurface bed 
material also tend to increase in size. 

Recognizing that armoring is a process that commonly 
leads to reduced scour in coarse-bed streams, a previously 
published concept relating to armoring was used to further 
evaluate variations in particle size between surface and 
shallow-subsurface bed material. Lisle and Madej (1992, 

p. 277) proposed the ratio of surface-to-subsurface bed- 
material particle sizes as a measure of armoring. The variation 
in bed-material size over a selected range of the particle-size 
distribution curve (D,, to D,,) was, therefore, computed for the 
surface and shallow-subsurface bed material associated with 
the pier-scour data collected for this study using the equation: 


(23) 
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where 
R is the bed-material size ratio of the 
corresponding surface and shallow- 
subsurface bed-material particle size; 

D, is the surface bed-material particle size, D, 
where i percent of particles are finer by 
weight, in millimeters or feet; and 

D. is the shallow-subsurface bed-material 
particle size, D, where i percent of particles 
are finer by weight, in millimeters or feet. 


Bed-material size ratios (median values for R,) for the 
coarser portion of the particle-size distribution curve (D,, to 
D,,) did not vary substantially, and ranged from about 1.1 to 
1.4 (fig. 21). The interquartile ranges (IR) also did not vary 
substantially and were all less than 1. Bed-material size ratios 
(median values for R,) for the finer portion of the particle-size 
distribution curve (D,, to D,,) were more variable than for 
bed-material size ratios for the coarser portion of the curve, 
and ranged from about 1.8 to about 4.0. The interquartile 
ranges also were more variable and ranged from about 1.2 to 
about 4.6. Thus, the greatest variability between surface and 
shallow-subsurface particle sizes, as defined by R., is at the 
fine-grained end of the particle-size distribution curve. 

The streambed upstream from the bridge likely was 
armored at most measurement sites because the surface bed 
material was more coarse when compared to the shallow- 
subsurface bed material, and because of predominant 
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Figure 21. Comparison of ratios of surface to shallow-subsurface bed-material particle size based on selected percentiles (D,) from 
the particle-size distribution curves for each site where pier-scour measurements were collected for this study. 


clear-water conditions (V < V,.,) and gradation coefficients 
greater than 1.4 (table 9). At the piers, the relative scour 
velocities (fig. 8) were great enough (V,/V"_., > 1.00) that the 
armor layers in the scour holes were mobilized and all sites 
were undergoing pier scour. Based on observed field condi- 
tions, and the relative coarseness of the bed-material particle 
sizes for the two layers (figs. 20 and 21), D,, and D,. particle- 
size data for the shallow-subsurface bed material underlying 
the armor layer (table 9) also were used to compute pier scour 
by the HEC-18-K,Mu equation (equation 5) with K, deter- 
mined by equation 11. 

Scour depths were computed using the D., and D,. shal- 
low-subsurface particle-size data for 94 of the 103 pier-scour 
measurements made for this study. Nine of the measurements 
either did not meet the criteria to be considered as a valid 
real-time scour measurement, or anomalous results were pro- 
duced when the D,, and D,. shallow-subsurface bed-material 
particle sizes were used in equation 11 to compute K,. Scour 
depths computed based on shallow-subsurface particle sizes 
were then compared to both measured scour depths (table 2) 
and scour depths computed based on the surface bed-material 
particle-size data (D,, and D,., table 9). Differences in D,, and 
D,, particle size between the surface and shallow-subsurface 
bed material did not cause substantial differences in computed 
scour (table 11). The frequency and magnitude of the residual 
pier-scour depth (fig. 22) associated with underpredicted and 
overpredicted scour was about the same regardless of which 
particle-size data were used with the HEC-18-K,Mu equa- 
tion. Computed scour was underpredicted in only 5 out of 
94 instances when surface bed-material particle-size data were 
used with the HEC-18-K,Mu equation and in only 3 out of 
94 instances when shallow-subsurface bed-material particle- 
size data were used. Positive (underpredicted) and negative 
(overpredicted) residuals of scour are noted as separate statis- 
tics in table 11 because of the implication each type of residual 
separately has to bridge foundation design. The median under- 
predicted (positive) residual scour using surface and shallow- 
subsurface data with the HEC-18-K,Mu equation (table 11) 
was 0.34 and 0.29 ft, respectively, and the maximum under- 
predicted residual scour was 0.52 and 0.34 ft, respectively. The 
median overpredicted (negative) residual scour using surface 
and shallow-subsurface data with the HEC-18-K,Mu equation 
(table 11) was -0.91 and -1.02 ft, respectively, and the greatest 
overpredicted residual scour was -4.41 ft using either bed- 
material layer. 

Computed scour generally was the same or slightly 
greater when shallow-subsurface bed-material particle-size 
data were used with the HEC-18-K,Mu equation because 
shallow-subsurface bed-material particle sizes (D,, and D,.) 
typically were slightly smaller than corresponding surface bed- 
material particle sizes. In addition, other factors contributed 
to the close agreement between pier-scour computations using 
the surface and shallow-subsurface bed-material size data. For 
example, of the 94 measurements made for this study satisfy- 
ing criteria for real-time scour at the pier (V, /V’.., > 1.00), 37 
had computed values of K, less than 0.4, regardless of whether 
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surface or shallow-subsurface bed-material data were used 
with equation 11. The minimum allowed value of K, (0.4) was 
substituted for the computed K, in both instances according to 
HEC-18 criteria (Richardson and Davis, 2001). This substi- 
tution resulted in identical computed scour depths for these 
measurements, regardless of which bed-material data were 
used. Furthermore, 34 of 94 measurements had both computed 
values of K, exceeding 0.4 for surface and shallow-subsurface 
bed-material data—the median difference in computed K, 
was 0.03. In addition, 21 of 94 measurements had at least one 
computed value of K, exceeding 0.4 when either surface or 
shallow-subsurface bed-material data were used with equa- 
tion 11—the median difference in computed K, was 0.05 after 
substituting 0.4 for the calculated K,, that was less than 0.4. 
Finally, substantial differences for the computed Mueller 
velocity-intensity ratio (V,,,,, equation 10) for some sites typi- 
cally resulted in much smaller differences between the result- 
ing values of K, computed by equation 11 using corresponding 
surface and shallow-subsurface data. The modest response 
of K, (dependent variable) to more widely varying values of 
Vx, (independent variable) is partly due to the small exponent 
(0.15) in equation 11, which indicates a relatively flat slope 
(weak relation) for the regression line relating K, to V,,,. In 
summary, despite vertical variations in coarse bed-material 
characteristics between the surface and shallow-subsurface 
bed-material layers, pier-scour depths computed using either 
set of characteristics were quite similar. The HEC-18-K,Mu 
equation was not very sensitive to vertical variations in bed- 
material characteristics within the active layer at coarse-bed 
sites where pier-scour data were collected for this study. 


Combined Effects of Surface and Shallow- 
Subsurface Bed-Material Characteristics 


In addition to examining the effects that differences 
in bed-material size and gradation between the surface and 
shallow-subsurface layers might have on predicted pier-scour 
depth, the combined effect that the bed-material characteristics 
of both layers might have on measured pier scour also was 
examined. The combined effect is of interest because if the 
critical velocity to initiate scour of bed material at the pier is 
exceeded, the surface and shallow-subsurface layers could 
interact and exchange bed material within the active layer until 
an immobile armor layer re-forms. Lagasse and others (1991) 
describe the armoring process as the interaction between the 
surface and shallow-subsurface layers of bed material. For 
gravel-bed streams where the bed is not composed of immo- 
bile armor but is paved so that coarse particles are frequently 
mobilized, Parker and others (1982, p. 546) indicate that the 
particle-size distribution of the bedload typically is more 
similar to the size distribution of the shallow-subsurface layer 
(subpavement) than to that of the surface layer (pavement). It 
is therefore conceivable that the bed-material size and grada- 
tion characteristics of the two layers in combination can affect 
pier scour. For example, substantial pier-scour depth might 
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Residual pier-scour depth is measured pier-scour 
depth minus computed pier-scour depth. 


Figure 22. Comparison of residual pier-scour depth based on pier-scour measurements and bed-material data collected for this study 
with pier-scour depth computed by the HEC-18-K Mu equation (Richardson and Davis, 2001). 


Table 12. Groups of pier-scour measurements based on classification of particle sizes of the 
surface and shallow-subsurface bed material. 


Bed-material size classification modified Number Group 

Combination from Lane and others (1947) of number 

Surface Shallow-subsurface measurements (fig. 23) 
1 Small cobble Small cobble 17 1 
2 Small cobble Very coarse gravel 9 2 
3) Very coarse gravel Small cobble 2 
4 Very coarse gravel Very coarse gravel 23 3 
5) Very coarse gravel Coarse gravel 14 4 
6 Coarse gravel Very coarse gravel 4 4 
7 Very coarse gravel Medium gravel 4 
8 Coarse gravel Coarse gravel 19 5 
9 Coarse gravel Medium gravel 1 5 
10 Medium gravel Medium gravel 7 Pe) 
11 Coarse gravel Very fine gravel 3 5 
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be expected where a thin layer or pavement of moderately 
coarse surface bed material like gravel is underlain by much 
finer subsurface bed material like sand (American Society of 
Civil Engineers, 1975, reprinted 1977, p. 182). In contrast, 
less scour might be expected when the bed material of the two 
layers includes coarse gravel and cobbles and the bed material 
of one layer has a greater degree of gradation than the other 
(fig. 4). 

The effects that differences in particle size between 
surface and shallow-subsurface bed-material layers might 
have on pier scour were investigated first. Examination of the 
differences in particle sizes between the two layers began by 
grouping the bed-material particle-size data for each layer 
into bed-material size classifications (table 3) based on the D,, 
particle size. The classifications ranged from small cobble to 
very fine gravel. The different combinations of surface and 
shallow-subsurface particle-size classifications for individual 
scour measurements were examined next. The 103 scour 
measurements collected for this study exhibited 11 combina- 
tions of bed-material size classification involving the surface 
and shallow-subsurface layers, and the measurements were 


initially grouped based on these 11 combinations (table 12). 
Many of the 11 combinations were characteristic of only a 
small number (four or less) of the 103 pier-scour measure- 
ments. Therefore, one or more combinations with only slightly 
different size characteristics were lumped together into five 
groups to allow development of envelope curves. The relation 
between relative velocity, based on the D., of the surface bed- 
material layer, and relative pier scour was examined for pier- 
scour measurements in each group (fig. 23). Based on the five 
envelope curves drawn to define the upper bounds of the data 
for each group, relative pier scour generally was less for mea- 
surements with coarser bed material like very coarse gravel or 
cobbles in the two layers (fig. 23). Few other conclusions are 
apparent from a comparison of the envelope curves. 

The relation between relative velocity, based on the aver- 
age D,, for the two bed-material layers, and relative pier scour 
also was examined and produced similar results. This similar- 
ity of results may be because surface and shallow-subsurface 
bed-material size classifications for both layers were the same 
or similar for both layers for most measurements (for example, 
small cobbles in both layers), even though surface bed- 
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(table 12) Surface layer Shallow-subsurface layer n is the number of measurements 
wana 1 Small cobble Small cobble y', is pier-scour depth from field measurements, in feet 
—i—. 2 Small cobble to very Small cobble to very 
coarse gravel coarse gravel b is effective width of the pier measured perpendicular to 
— 3 Very coarse gravel Very coarse gravel the approach flow, in feet 
ie, as 4 Very coarse gravel Very coarse gravel to V, isthe measured approach velocity directly upstream 
to coarse gravel medium gravel from the pier, in feet per second 
—_—e— 5 Coarse gravel to Coarse gravel to very 


medium gravel 


fine gravel 


is the critical velocity for incipient motion of surface bed 
material upstream from the bridge for the D., particle size, 
in feet per second 


material size typically was larger than the shallow-subsurface 
bed-material size. There were no sites where the difference in 
bed-material particle size between the surface or armor layer 
and shallow subsurface was substantial (for example, small 
cobbles underlain by very fine gravel or sand); thus, more 
detailed assessment of such widely contrasting bed-material 
compositions was not possible. 

Differences in gradation of bed material in the sur- 
face and shallow-subsurface layers and the relation of these 
differences to pier scour also were investigated using two 
approaches. First, the differences in relations between grada- 
tion coefficients for each layer and relative pier scour were 
compared by examining envelope curves for each set of mea- 
surements obtained for this study (fig. 24). Maximum relative 
pier scour corresponded to a gradation coefficient of about 
2.3 for the surface bed material and about 4.7 for the shallow- 
subsurface bed material (fig. 24). This fairly narrow range in 
gradation coefficients associated with maximum pier scour 
for the two layers is similar to other studies that showed the 
narrow range in gradation coefficients associated with maxi- 
mum pier scour based only on surface bed material (Ettema, 
1980; Baker, 1986). A composite envelope curve (fig. 24) was 
generated as a means to account for the variability in 0, as it 
pertains to the computation of pier-scour depth when results of 
multiple pebble counts or sieve analyses involving both layers 
are combined. The composite gradation coefficient (9, ) is the 
geometric mean of the corresponding surface and subsurface 
o, values. 
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For the second approach, the variation in gradation over 
the range of the particle-size distribution curve was defined 
using the following relation: 


(24) 


where 
R, is the relative gradation or the ratio of 
shallow-subsurface-to-surface bed-material 
gradation coefficient, dimensionless; 

c.. is the gradation coefficient of the shallow- 
subsurface bed material, dimensionless; 
and 

c,, is the gradation coefficient of the surface bed 
material, dimensionless. 


The relative gradation (R,) was calculated for each of the 

103 pier-scour measurements collected for this study and 

the measurements were ranked by R, from largest to small- 
est. Selected percentiles (99th, 75th, 50th, 25th, and 10th) 
were then determined and the relative pier scour and relative 
velocity for each measurement in the percentile were plotted 
to examine the relation between relative pier scour and relative 
velocity for each percentile. Envelope curves bounding the 
measurements in each percentile were drawn. Each of the five 
envelope curves (fig. 25), therefore, bounds a group of pier- 
scour measurements with RX, values that equal or exceed a par- 
ticular value (smallest value) in each percentile. For example, 
10 percent of the measurements had R, values that equaled or 


EXPLANATION 
Envelope curves and data points 
—e— Surface gradation coefficient, G,. 
——*—.__ Shallow-subsurface gradation coefficient, Se 
---@--- Composite gradation coefficient, G,. based on 
geometric mean of G,. and o.., 
Variables 
n is the number of measurements 
y', is pier-scour depth from field measurements, 
in feet 
b is effective width of the pier measured perpendicular 
to the approach flow, in feet 
G,, is the surface bed-material gradation coefficient, 
: dimensionless 
ae is the shallow-subsurface bed-material gradation 
. coefficient, dimensionless 
OF, is the composite gradation coefficient or the geometric 
mean of the corresponding surface and shallow- 
subsurface G, values, dimensionless 
oO is the bed-material gradation coefficient, dimensionless 


Figure 24. Comparison of relative pier scour to gradation coefficient for surface and shallow-subsurface bed-material data collected 


for this study. 
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exceeded 3.18, so 10 data points (R, = 3.18, fig. 25) were used 
to define the 10th-percentile envelope curve. Maximum rela- 
tive scour associated with each envelope curve is concentrated 
at about the same relative velocity (V/V... ~ 0.75 to 0.8). The 
five envelope curves (fig. 25) appear to converge when rela- 
tive velocities are high enough for live-bed scour conditions 
to be reached (V/V. > 1.00). Based on these curves, when 
the difference between the shallow-subsurface and surface 
bed-material gradation coefficients increases (higher values of 
R,), pier scour decreases. Earlier studies reported by Mueller 
(1996) and Mueller and Wagner (2005) indicated that when 
the gradation coefficient of the surface bed material increased, 
pier scour was reduced. 

Higher variability in gradation coefficients would seem to 
indicate the presence of a higher fraction of fine-grained par- 
ticles with expected deeper scour before armoring. Increased 
scour associated with finer bed material may, however, be 
more than offset if sufficiently large particles also are avail- 
able to hide the smaller particles (Einstein, 1950; Parker and 
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Klingeman, 1982), limiting overall mobility of bed material 
and allowing more rapid armoring. 

Finally, the combined influence of relative differences in 
surface and shallow-subsurface bed-material size and grada- 
tion on pier scour was investigated by defining a bed-material 
variability index (/)): 


D D 
i. = 50,5 ny 95,5 £58 (25) 
D5 5s Ds 5s Ons 
Substitution based on equations 23 and 24 yields: 
1, = Rag t Rye +k, (26) 


where all terms are as previously defined. The three variables 
on the right hand side of equation 26 were derived from 
variables (D,,, D,., and 9,) shown to be related to pier scour in 
coarse-bed streams as noted in this section and other studies 
cited in the section “Development of Pier-Scour Equations for 


4 EXPLANATION 
Envelope curves and data points 


99th percentile; Rg greater than or equal to 0.96 


) S=— 75th percentile; Rg greater than or equal to 1.24 
+ —eo— 50th percentile; Rg greater than or equal to 1.95 
7, —_ 25th percentile; Rg greater than or equal to 2.46 
] —— 10th percentile; Rg greater than or equal to 3.18 


Some data points associated with an envelope 
curve have coordinates identical to data points 

4 associated with other envelope curves. In those 
Al instances, the data points were slightly offset from 
each other in the horizontal, vertical, or both 

4 directions for clarity. 


zl Variables 
| n is the number of measurements 
i Vo is the approach velocity just upstream from the 
q pier, in feet per second 
| Vos0 is the critical velocity for incipient motion of 
| surface bed material upstream from the 
bridge, in feet per second 
1.3 1.4 P ; ; ; 

Rg is the relative gradation or the ratio of shallow 
subsurface-to-surface bed-material 
gradation coefficients, dimensionless 

Sgs is the surface bed-material gradation coefficient, 
dimensionless 

Sg ss is the shallow-subsurface bed-material 


gradation coefficient, dimensionless 


Figure 25. Relation between relative pier scour and relative velocity based on surface and shallow-subsurface bed-material gradation 


differences for pier-scour measurements collected for this study. 


Coarse-Bed Streams.” These variables have about the same 
magnitude and, therefore, have a similar influence on the val- 
ues of /. Based on research cited in the section “Development 
of Pier-Scour Equations in Coarse-Bed Streams,” that showed 
pier scour to be less for coarse-size bed material with higher 
gradation coefficients, perfect homogeneity (where R.,, R,., 
and R, each equal 1.00 in equation 26) between surface and 
shallow-subsurface bed-material size and gradation (/, = 3.00) 
might be expected to result in greater relative pier scour than 
a condition where, for example, the relative homogeneity of 
size and gradation varied by a factor of two for each ratio 
(I, = 6.00). 

The bed-material variability index (/,) was calculated 
for each of the 103 measurements collected for this study and 
the measurements were ranked by /, from largest to smallest. 
Selected percentiles (99th, 75th, 50th, 25th, and 10th) were 
determined and the relative pier scour and relative velocity for 
each measurement in the percentile were plotted to examine 
the relation between relative pier scour and relative velocity 
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Figure 26. Relation between relative pier scour and relative velocity based on 
surface and shallow-subsurface bed-material size and gradation differences for 


pier-scour measurements collected for this study. 
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for each percentile. Envelope curves bounding the measure- 
ments in each percentile were drawn. Each of the five enve- 
lope curves (fig. 26), therefore, bounds a group of pier-scour 
measurements with J, values that equal or exceed a particular 
value (smallest value) in each percentile group. For example, 
10 percent of the measurements had /, values that equaled or 
exceeded 6.48, so 10 data points (/, > 6.48, fig. 26) were used 
to define the 10th-percentile envelope curve. Maximum rela- 
tive scour associated with each envelope curve is concentrated 
at about the same relative velocity (V/V _., ~ 0.80). The five 
envelope curves (fig. 26) appear to converge when relative 
velocities are high enough for live-bed scour conditions to be 
reached (V,/V_., > 1.00). Based on these curves, as differences 
in particle size and gradation between the surface and shallow- 
subsurface bed-material layers increase (higher values of J), 
pier scour decreases. These curves also exhibit the potential to 
refine estimates of maximum relative pier scour using the bed- 
material variability index or a similar variable. 
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ile; J, greater than or equal to 4.60 
25th percentile; /, greater than or equal to 5.52 
10th percentile; 7, greater than or equal to 6.48 
Some data points associated with an envelope curve 
have coordinates identical to data points associated 
with other envelope curves. In those instances, the 


data points were slightly offset from each other in the 
horizontal, vertical, or both directions for clarity. 


Variables 


is the number of measurements 


is pier-scour depth from field measurements, 
in feet 


is effective width of the pier measured perpendicular 
to the approach flow, in feet 


is the bed-material variability index, dimensionless 


13 14 Vi is approach velocity directly upstream from the pier, 
in feet per second 
V so is the critical velocity for incipient motion of the 
: surface bed material upstream from the bridge 
for the Dsq particle size, in feet per second 
Dis) is the surface bed-material particle size, D, where 
i percent of particles are finer by weight, in 
millimeters or feet 
Dass) is the shallow-subsurface bed-material particle size, 
D, where i percent of particles are finer by 
weight, in millimeters or feet 
Sg is the surface bed-material gradation coefficient, 
dimensionless 
Sg ss is the shallow-subsurface bed-material gradation 


coefficient, dimensionless 
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Limitations of the Study 


Field measurement of scour at bridge piers requires the 
use of specialized equipment and methods to obtain accurate 
data. The collection of real-time scour data also requires that 
measurements be made during high-flow conditions. Although 
methods for conducting scour-depth measurements gener- 
ally are described in various USGS publications, innovation, 
physical effort, and judgment are required in obtaining usable 
pier-scour data under difficult field conditions. A series of 
measurement and interpretive steps are needed to compute the 
measured pier-scour depth. Except for three sites where con- 
tinuous scour measurements were made over the entire runoff 
period by using fixed instrumentation, most measurements for 
this study were conducted over a few hours on a single day. 
Continuous transducer soundings made at one site indicated 
that maximum pier scour was attained several times during a 
2-month period in 2006 for flows that ranged from less than 
to about bankfull discharge, which provide some insight on 
the site-specific effects of high-flow duration on pier scour. 
High-flow duration effects on pier scour, however, were not 
addressed beyond this one example. 

Runoff conditions during the study mostly were lim- 
ited to bankfull discharge approximating the 1.5- to 2-year 
recurrence-interval flood. Lack of higher flows having greater 
velocity intensities may indicate that scour depths were limited 
by hydrologic conditions. Still, scour measurements made 
mostly under clear-water conditions associated with maximum 
scour involved relatively high velocity intensities. In fact, 
measured velocities at one-half of study sites were at least 
70 percent of the critical velocity for incipient motion of the 
bed material upstream from the bridge, which might indicate 
that scour was measured near the threshold between clear- 
water and live-bed scour conditions. Laboratory research using 
fine-grained material and steady-flow conditions indicates that 
this threshold is where maximum pier scour likely occurs. 
Furthermore, historical pier-scour data from Montana with 
peak streamflows of greater recurrence interval (up to about 
100 years) plotted below envelope curves (figs. 10, 11, and 14) 
that bound the data collected for this study and BSDMS data. 
Pier-scour depths measured in this study also may have been 
limited by streambed armoring and unsteady flow conditions. 

It is critical for bed-material samples to be representative 
of the actual streambed materials if these data are to be related 
to maximum pier-scour depth and used to improve equations 
to predict pier scour. A single bed-material sample considered 
to be representative of overall conditions was obtained for 
most sites. At 10 percent of the sites, multiple sets of pebble 
counts (3 to 4, table 4) were obtained. Results indicated that 
the standard deviation of the particle-size was small in most 
instances when compared to the corresponding mean particle 
size at each percentile analyzed. 

Corresponding data for surface and shallow-subsurface 
bed-material characteristics obtained at all sites provided 
an indication of the vertical variability of the particle size. 
Relations between relative pier scour and relative velocity 


involving surface and shallow-subsurface bed-material char- 
acteristics were limited to the data obtained in this study; no 
comparable data were found in the BSDMS database. Thus, 
surface-to-shallow-subsurface relations need to be considered 
as preliminary. Additional data from bridge sites with a wider 
range in bed-material size and gradation are needed to verify 
the relations presented and perhaps use this information to 
improve scour-prediction equations for use at sites with simi- 
lar bed-material characteristics. An independent dataset also 
would be needed to verify the validity of surface-to-shallow- 
subsurface relations proposed in this study. 

Although differences between surface and shallow- 
subsurface bed-material size and gradation resulted in no 
substantial difference in computed scour depths, the spatial 
variability in bed-material composition needs to be considered 
in assessing the site-specific effect that bed-material size and 
gradation differences have on predicted pier scour. Further- 
more, there may be combinations not investigated in this study 
where the bed-material differences in the two layers are more 
substantial and lead to greater differences in predicted scour 
depth. Future investigations of relations between bed-material 
characteristics and pier scour will require detailed descrip- 
tion of sampling location (areally, vertically, or both) for the 
bed-material data collected. Some difficulty is likely to be 
encountered in collecting subsurface bed-material samples at 
sites where the streambed is always under a substantial depth 
of water or velocities are excessive, which is common in some 
regions of the United States. Conclusions presented for this 
study are generally limited to the range of hydraulic conditions 
and bed-material characteristics demonstrated in the data col- 
lected for this study and BSDMS data. 


Summary and Conclusions 


Determination of pier-scour potential is an important 
consideration in the hydraulic analysis and design of highway 
bridges that cross streams, rivers, and other waterways in the 
United States. A primary goal of ongoing research in the field 
of bridge scour is to improve equations currently used (2011) 
so that scour depth is neither underpredicted nor excessively 
overpredicted. Research efforts, therefore, continue to evaluate 
and improve the accuracy of pier-scour equations. Since the 
early 1990s, research has expanded to include field measure- 
ments of pier scour during high-flow conditions. Based on 
these measurements, scour depth for piers in noncohesive, 
nonuniform streambeds with a mixture of sand, gravel, 
cobbles, and boulders (coarse-bed streams) generally is less 
than the scour depth in finer-grained (mostly sand) streambeds 
under similar hydraulic conditions. The difference in scour 
depth is attributed to formation of an armor layer. This finding, 
supported by field measurements collected by the U.S. Geo- 
logical Survey (USGS), has led to modifications in the scour 
equations recommended by the Federal Highway Administra- 
tion (FHWA). Pier-scour data collected by the USGS were 


used to develop a bed-material correction factor called K,. 
Collection of additional coarse-bed pier-scour data was recom- 
mended to further validate the method for computing K, and 
to provide additional data for continued evaluation of pier- 
scour equations. The Montana Department of Transportation 
(MDT) was interested in pier-scour prediction in coarse-bed 
streams because coarse-bed streams are common in Montana. 
Furthermore, the USGS has collected bridge-scour data and 
conducted bridge-scour investigations since the early 1990s 
in cooperation with MDT. Consequently, the USGS and MDT 
began a cooperative study in 2001 to investigate pier scour in 
coarse-bed streams in Montana. 

This report presents the results of a study of pier scour in 
coarse-bed streams during 2001-07 in the mountain and foot- 
hill regions of three major river basins in Montana: the Mis- 
souri, Yellowstone, and Clark Fork of the Columbia. Drainage 
areas for the streams at bridge sites where measurements were 
collected ranged from about 3 mi? to almost 20,000 mi’. The 
report also describes how coarse-bed pier-scour measurements 
were collected, shows the extent that the coarse portion of the 
national pier-scour database was expanded, discusses how 
these new data were used to evaluate the relative accuracy of 
various equations for predicting scour in coarse-bed streams, 
and demonstrates how differences in size and gradation 
between surface bed material and bed material underlying the 
surface layer (shallow-subsurface bed material) might relate to 
pier scour. 

The equation widely adopted for conducting pier-scour 
analysis and design in the United States was developed at 
Colorado State University in the mid—1970s. A modified ver- 
sion of the equation was later published by FHWA in Hydrau- 
lic Engineering Circular No. 18 (HEC-18) and is called the 
HEC-18 equation. The current (2011) version of the HEC-18 
pier-scour equation recommended by the FHWA includes 
the computation of K,, and is referred to in this report as the 
HEC-18-K,Mu equation. Pier-scour research has continued 
and new equations have been proposed that include the K, 
correction factor. 

Methods based on limited-detail field procedures 
described in literature and equipment ranging from simple 
scour rods to complex hydroacoustic instrumentation were 
used to obtain pier-scour measurements. Data collected in 
the field for this study included information needed to mea- 
sure pier scour along with bed material and ancillary data 
needed to compute scour using equations. Bed-material sizes 
used in scour-prediction equations were determined from the 
surface bed material that generally reflected the armor layer. 
The potential interaction between the surface and shallow- 
subsurface bed-material layers with respect to pier scour 
and armoring also was of interest for this study. Therefore, 
shallow-subsurface bed material also was sampled and data 
were compared to surface bed-material data for each site. Sur- 
face and shallow-subsurface bed-material measurements were 
obtained by pebble counts and sieve analyses. Bed-material 
data were used to construct particle-size distribution curves, 
which were used to conduct various coarse-bed pier-scour 
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analyses. Surface and shallow-subsurface bed material 
sampled for this study was used to evaluate the sensitivity of 
computed scour to size and gradation differences between the 
two layers, and these data for each layer were used to assess 
the effects of bed-material characteristics on pier scour. 

Pier-scour measurements were made during 2001 through 
2007 at 59 bridge sites on coarse-bed streams. Of those mea- 
surements, 103 satisfied criteria requiring that the measured 
approach velocity (V,) exceed the critical velocity (V’..,) to 
initiate scour of bed material at the piers and were analyzed 
for this study. Incipient motion computations based on the D,, 
particle size indicated that 96 out of 103 pier-scour measure- 
ments were made under clear-water scour conditions, when 
the streambed upstream from the bridge is stable and there is 
no substantial incoming sediment supply to the bridge open- 
ing. Of the measurements made, 50 percent had an approach 
velocity that equaled or exceeded 70 percent of the critical 
velocity (V,,,) for incipient motion of bed material upstream 
from the bridge. This might indicate that scour was measured 
near the threshold between clear-water and live-bed scour 
(V,/V.., = 1.00) where maximum pier scour was shown in 
laboratory studies. Despite reasonably high approach veloci- 
ties, measured pier scour was less than a couple of feet in 
depth for 70 percent of the measurements made at the coarse- 
bed sites. Particle size of sampled surface bed material was 
predominantly very coarse gravel based on the D., particle 
size. Surface bed material at one-half of the sites lacked sub- 
stantial gradation, as indicated by gradation coefficients (0,) 
less than about 2.0. About 60 percent of the sites were located 
on high-gradient streams based on water-surface slopes that 
exceeded 0.002. The approach depth of flow (y,) in relation 
to the effective pier width (b) was not substantial, with flow 
depth typically less than twice the effective pier width. The 
effects of high-flow duration on pier scour were not investi- 
gated because most scour data for a site were collected in a 
single day. Continuous scour measurements made at three sites 
using fixed acoustic transducers, however, indicated that scour 
typically increased as streamflow increased. 

Pier-scour data collected for this study were compared 
to pier-scour data from other USGS studies that were used 
to develop pier-scour equations for coarse-bed streams. 
Historical data were retrieved from the Bridge Scour Data 
Management System (BSDMS). The 508 measurements in 
the BSDMS database collected between 1965 and 1998 were 
initially screened down to 272 measurements (BSDMS All) 
with bed material that was noncohesive and composed of sand 
and coarser particles, and then to 90 measurements (BSDMS 
Coarse). The BSDMS Coarse dataset is comparable with data 
collected for this study based on FHWA criteria for D,, and 
D,, particle size that defines coarse bed material. The median 
D,, and D,, bed-material size for data collected for this study 
are about 1.8 and 1.6 times (respectively) the size of D,, 
and D,. in the BSDMS Coarse dataset. The maximum value 
of relative pier scour (y’ /b) in the study was 1.00, and the 
maximum values observed in the two BSDMS datasets (All, 
Coarse) were 2.09 and 1.25, respectively. Data from this study 
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also were compared to historical pier-scour measurements 
made in Montana by USGS during large floods. Nine pier- 
scour measurements made at three sites were for peak stream- 
flows approximating the 100-, 10-, and 5-year recurrence 
intervals, which are much greater than any recurrence intervals 
associated with measurements made for this study and provide 
an independent check on the envelope curves bounding the 
upper end of both the data collected for this study and BSDMS 
data combined. 

Pier-scour depth was computed for the measurements 
collected for this study using the HEC-18 equation with no K, 
correction factor and five pier-scour equations that use K, for 
the armoring effect of coarse bed material. Computed scour 
depths were compared to measured scour depths using scatter 
plots and boxplots of residuals to evaluate the reliability of the 
equations to predict pier scour. The HEC-18-K,Mu equation 
was the best equation for predicting pier-scour depth in coarse- 
bed streams because the equation generally predicted pier 
scour in closer agreement to measured scour than the other 
equations used for computing pier-scour depth. Furthermore, 
the magnitude of the residuals of underpredicted and overpre- 
dicted scour depth was the lowest using the HEC-18-K,Mu 
equation. 

The effect of bed-material size and gradation on scour 
depth was investigated using data collected for this study, 
the BSDMS All dataset, and the BSDMS Coarse dataset to 
confirm and extend the findings of previous studies. Previ- 
ous research showed that the equation for computing K, can 
produce anomalous results (negative velocity ratio) for rare 
combinations of D,, and D,.. None of the surface bed-material 
particle-size data collected for this study produced anomalous 
results in computing K,. Regression analysis conducted in 
this study showed that the D,, and D,, particle sizes have a 
positive relation with each other that is statistically significant 
(r? = 0.73 and p < 0.001). In other words, for the data col- 
lected for this study, sites that have a larger D,, particle size 
also tend to have a larger D,. particle size. This observation is 
similar to the observation of previous research about the bed- 
material data originally used to develop the equation for the K, 
correction factor. 

The relation between relative pier scour (y’,/b) and rela- 
tive velocity (V, /V_.,) defined by an envelope curve for the 
data collected for this study and the BSDMS Coarse dataset 
were compared to relations observed in laboratory data. The 
envelope curve for data used for this study displays an earlier 
peak, with a decrease in relative pier scour followed by an 
asymptotic rise; the first peak in the curve occurs at a lower 
relative velocity than the clear-water and live-bed threshold 
(V,/V.,, = 1.00). The magnitude of relative pier scour (maxi- 
mum of 1.25) is much less than the published maximum value 
of 2.4 obtained from laboratory data. Furthermore, equilib- 
rium live-bed scour is about 18 percent shallower than the 
maximum pier scour under clear-water conditions, a larger 
reduction than the 10 percent indicated in literature. Armoring 
associated with nonuniform coarse-bed material and unsteady- 
flow conditions in the field probably accounts for differences 


between envelope curves developed for this study and those 
from previous studies conducted in the laboratory using fine- 
grained material and steady-flow conditions. 

Previous studies showed that surface bed-material grada- 
tion affected pier-scour depth, therefore, the relation between 
gradation coefficient (o,) and relative pier scour was inves- 
tigated using the pier-scour measurements collected for this 
study. Gradation coefficients associated with surface bed mate- 
rial collected for this study ranged from 1.48 to 4.14, with a 
median of 2.01. Bed-material data collected for this study had 
sufficient gradation so that streambed armoring was possible 
according to criteria developed in previous research. As shown 
in previous research, for data collected for this study and 
the BSDMS Coarse dataset, maximum relative pier scour is 
greater when o, is less than or equal to 2.5. Pier-scour mea- 
surements collected for this study and BSDMS confirm a gen- 
eral lack of data for coarse-bed sites with o, greater than about 
2.5. The narrow range in gradation coefficients displayed by 
all available data and the nonuniform conditions present in 
the field continue to complicate efforts to isolate the effects of 
gradation on pier scour. 

Since the mid—1960s, pier-scour research has focused 
primarily on surface bed material effects on pier scour. 
Surface and shallow-subsurface bed material were collected 
and analyzed for both layers during this study. Median D,, 
particle sizes for 103 surface and shallow-subsurface bed- 
material samples were about 49 mm and 32 mm, respectively, 
and median D,. particle sizes for 103 surface and shallow- 
subsurface bed material samples were 122 mm and 86 mm, 
respectively. The median gradation coefficients for surface 
and shallow-subsurface bed material associated with measure- 
ments made for this study were 2.01 and 4.14, respectively. 
The combined effect of the larger surface bed-material size 
relative to shallow-subsurface bed-material size and the 
wider range of gradation coefficients in shallow-subsurface 
bed material is typical for the pier-scour measurements made 
for this study and may contribute to reduced pier scour in 
these coarse-bed streams. A significant relation was indicated 
between surface and shallow-subsurface bed-material size, 
for the D,, and D,, particle sizes (for D.,, 7° = 0.690 and 
p < 0.001; for D,., 7° = 0.647 and p < 0.001). In other words, 
as the particles in the surface bed material increase in size, the 
particles in the shallow-subsurface bed material also tend to 
increase in size. Based on observed field conditions, and the 
relative coarseness of the bed-material particle sizes for the 
two layers, particle-size data for the shallow-subsurface bed 
material underlying the armor layer also were used to compute 
pier scour by the HEC-18-K,Mu equation. Scour depths com- 
puted based on shallow-subsurface particle sizes were then 
compared to both measured scour depths and scour depths 
computed based on the surface bed-material particle-size data. 
Differences in D,, and D,, particle size between the surface 
and shallow-subsurface bed material did not cause substantial 
differences in computed scour. The frequency and magnitude 
of the residual pier-scour depth associated with underpredicted 
and overpredicted scour was about the same regardless of 


which particle-size data were used with the HEC-18-K ,Mu 
equation. 

In addition to examining the effects that differences 
in bed-material size and gradation between the surface and 
shallow-subsurface layers might have on predicted pier-scour 
depth, the combined effect that the bed-material characteristics 
of both layers might have on measured pier scour also was 
examined. Examination of the differences in particle sizes 
between the two layers began by grouping the bed-material 
particle-size data for each layer into classifications that ranged 
from small cobble to very fine gravel. Different combinations 
of surface and shallow-subsurface bed-material size classifica- 
tions were used to group the pier-scour measurements, and the 
relation between relative pier scour and relative velocity was 
examined for each of these groups of measurements. Based on 
the five envelope curves drawn to define the upper bounds of 
the data for each group, relative pier scour was generally less 
for pier-scour measurements with coarser bed material like 
very coarse gravel or cobbles in the two layers. There were no 
sites where the difference in bed-material particle size between 
the surface or armor layer and shallow-subsurface layer was 
substantial (for example, small cobbles underlain by very fine 
gravel or sand); thus, more detailed assessment of such widely 
contrasting bed-material compositions was not possible. 

Differences in gradation of bed material in the sur- 
face and shallow-subsurface layers and the relation of these 
differences to pier scour also were investigated using two 
approaches. First, the differences in relations between grada- 
tion coefficients for each layer and relative pier scour were 
compared by examining envelope curves for each set of 
measurements obtained for this study. Maximum relative pier 
scour corresponded to a gradation coefficient of about 2.3 
for the surface bed material and about 4.7 for the shallow- 
subsurface bed material. A composite curve was proposed 
that incorporates the separate relation of both layers. For the 
second approach, the variation in gradation over the range of 
the particle-size distribution curve was defined using ratios of 
shallow-subsurface-to-surface gradation coefficient (rela- 
tive gradation, R_) for each of 103 pier-scour measurements 
collected for this study. The measurements were ranked by 
R, from largest to smallest, selected percentiles were deter- 
mined, and the relative pier scour and relative velocity for 
each measurement in the percentile were plotted to examine 
the relation between relative pier scour and relative velocity 
for each percentile. Envelope curves bounding the measure- 
ments in each percentile were drawn. Maximum relative scour 
associated with each envelope curve was concentrated at about 
the same relative velocity. Based on these curves, when the 
difference between the shallow-subsurface and surface bed- 
material gradation coefficients increases (higher values of R.), 
pier scour decreases. 

Finally, the combined influence of relative differences 
in surface and shallow-subsurface bed-material size and 
gradation on pier scour was investigated by defining a bed- 
material variability index (/,) based on ratios of particle sizes 
and gradation coefficients between layers. The bed-material 
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variability index was calculated for each of the 103 measure- 
ments collected for this study and the measurements were 
ranked by J, from largest to smallest. Selected percentiles were 
determined and the relative pier scour and relative velocity for 
each measurement in the percentile were plotted to examine 
the relation between relative pier scour and relative velocity 
for each percentile. Envelope curves bounding the measure- 
ments in each percentile were drawn. Maximum relative scour 
associated with each envelope curve was concentrated at about 
the same relative velocity. Based on these curves, as differ- 
ences in particle size and gradation between the surface and 
shallow-subsurface bed-material layers increase (higher values 
of I), pier scour decreases. 

Field measurement of scour at bridge piers requires the 
use of specialized equipment and methods to obtain accurate 
data under difficult field conditions, and interpretive steps are 
involved to compute the measured pier-scour depth. Runoff 
conditions during the study mostly were limited to bankfull 
discharge approximating the 1.5- to 2-year recurrence-interval 
flood, which have a 67- and 50-percent chance, respectively, 
of happening in any year. Lack of higher flows having greater 
velocity intensities may indicate that scour depths were lim- 
ited by hydrologic conditions. Pier-scour depths measured in 
this study also may have been limited by streambed armoring 
and unsteady flow conditions. A single bed-material sample 
considered to be representative of overall conditions was 
obtained for most sites. Relations between relative pier scour 
and relative velocity involving surface and shallow-subsurface 
bed-material characteristics were limited to the data obtained 
in this study; no comparable data were found in the BSDMS 
database. Thus, surface-to-shallow-subsurface relations need 
to be considered as preliminary. Conclusions presented for this 
study are generally limited to the range of hydraulic conditions 
and bed-material characteristics demonstrated in the data col- 
lected for this study and the BSDMS data. 
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